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L. RIASSUNTO 



Viene descritto un equalizzatore per riceviton di stazioni radio base di urr 
sistema radiomobile cellulare (GSM o DCS) che governa un array ad apertura 
sintetica per I'inseguimento dei mobili, operante in presenza di canali tempo-vananti 
e/o affetti da interferenti non-stazionari. Un insieme di array sistolici, posto a monte di 
un processore di Viterbi, minimizza I'errore quadratico medio che nasce dal confronto 
tra i simboli ricevuti, fiitrati spazialmente, ed i simboli della sequenza di training 
convoluti con i campioni della risposta impulsiva stimata del canale. La 
minimizzazione deirerrore sfrutta in mpdo nuovoia "Displacement Theory" applicata 
alia decomposizione in forma triangolare delle matrici ^ algebriche mediante 
aggiornamento ricorsivo dei cosiddetto "fattore di Cholesky". E prevista una variante 
per il caso di improwisi e forti interferenti, in cui viene utilizzato un rivelatore dei picchi 
di energia dei segnali sulle antenne per comandare raggiornamento anticipato deila 
sola componente spaziale di una matrice strutturata ottenuta mediante la suddetta 
decomposizione (fig. 12). 
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DESCRIZIONE 
CAMPO PI APPUCAZIONE Ml 5 B A 1 2 8 0 

La presente invenzione si riferisce al campo della telefonia radiomobile cellulare 
e piu precisamente ad un metodo di equalizzazione spaziale e temporale, a 
convergenza rapida, per la cancellazione di interferenti isofrequenziali stazionari, e 
non, in segnali ricevuti da antenne a schiera di stazioni radio base di un sistema 
radiomobile cellulare, e relativo equalizzatore ad array sistolici. In modo particolare, 
I'invenzione k riferita a sistemi radiomobiii che utilizzano una breve sequenza di 
training nei burst di trasmissione. 

Nei sistemi telefonici radiomobiii, la limitazione principale alia capacita di traffico 
offerta al pubblico servizio 6 costituita dalla scarsita di banda a disposizione. Sono 
quindi particolarmente apprezzate tutte quelle soluzioni tecnologiche che consentono 
di aumentare la suddetta capacity o efficienza spettrale, mantenendo una 
predeterminata qualita di servizio in un'area comunque estesa. Sono state alio scopo 
^ sviluppate opportune strategie di copertura cellulare, unitamente a tecniche di 

multiplazione in frequenza FDMA (Frequency Division Multiple Access), e/o nei tempo 
TDMA (Time Division Multiple Access), e piu recentemente nello spazio SDMA (Space 
Division Multiple Access) grazie airutilizzo di antenne a schiera cosiddette intelligenti 
abbinate ai ricevitori. 

La tecnica di copertura cellulare consiste nei suddividere I'area di servizio in 
zone, dette celle, ciascuna delle quali 6 servita da : una propria stazione radio base BS 
(Base Station). Ogni cella ha a disposizione un proprio gruppo di canali radio che pud 



essere riutilizzato in altre celle opportunamente distanziate, I'attenuazione dei segnah 
radio 6 infatti proporzionale alia distanza con un esponente, o coefficiente di 
propagazione, in genere maggiore di 3. Uinsieme complete delle celle che utilizzano 
tutti i canali radio a disposizione del sistema, senza alcun riuso degli stessi, e detto 
cluster di celle. Essendo I'efficienza spettrale inversamente proporzionale ali'area di 



cella, i sistemi microcellulari sono Tapproccio piu comune per aumentare I'efficienza 
spettrale in aree urbane o suburbane ad alta densita di traffico. Al diminuire dell'area 
di cella i problem! associati a continui hand-over diventano sempre piu critici, 
conseguentemente vengono in genere preterite strutture miste (microcelle 
sovrapposte a macrocelle ad ombrello) per gestire contemporaneamente stazioni 
mobili MS (Mobile Station) veloci in aree scarsamente popolate, e stazioni mobili lente 



in aree ad alta densita. 

Al fine di incrementare I'efficienza spettrale, un primo approccio alternativo alia 
riduzione della dimensione della cella, consiste nella riduzione del numero C di celle 
per cluster, dove C 6 detta anche dimensione del cluster, il che equivale ad 
aumentare la frazione dello spettro totale gestito dalla singola BS. Questa strategia 6 
nota come RCS (Reduced Cluster Size), il cui limite teorico e quello di rendere 
disponibile I'intero spettro a tutte le celle, reaiizzando cluster di una sola cella. 

Un secondo approccio consiste nell'allocare un numero maggiore di utenti nello 
spettro gi& a disposizione nella cella, allocando piu mobili sullo stesso canale, 
caratterizzato dalla stessa portante e dallo stesso time slot e, quindi, senza intervenire 
sulla dimensione del cluster. Questa seconda strategia 6 nota come SCFR (Same 
Cell Frequency Reuse) ed 6 basata sulla multiplazione a divisione di spazio SDMA 
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entro la cella. Nei casi pratici essa viene realizzata in due modi distinti, di cui un pTftflb 
consiste nella divisione della cella in settori piu piccoli mediante una disposizione di 
antenne con diagramma di radiazione fisso (in genere con apertura a 120°), ed un 
secondo, e piu moderno, consiste nel dotare le BS di rispettive antenne multifascio, o 



mobili con fasci direttivi relativamente stretti; cio implica che tali fasci debbano 



La formale differenza tra I'impiego di una tecnica SCFR e la suddivisione della 
cella stessa in settori piu piccoli, risiede nella capacita di allocare gli utenti nei rispettivi 
canali. Con la settorializzazione e possibile gestire piu utenti sullo stesso canale solo 
se la loro poslzione sul territorio rientra in different! settori fissi di copertura radio. 
Questo com porta, in uno scenario caratterizzato da utenti in movimento, la necessity 
di continui handover. Una struttura SCFR riesce invece ad allocare nella stessa cella 
piu utenti cocanale, senza dover partizionare la cella e ricorrere quindi a frequenti 
handover. Owiamente, anche in quest'ultimo caso la BS deve poter distinguere I'uno 
dairaltro gli utenti cocanale, in modo da formare in ricezione degli opportuni 
diagrammi di radiazione, ciascuno in grado di tracciare il rispettivo utile e reiettare gli 
interferenti isofrequenziali, corrispondenti ai rimanenti utenti cocanale 
simultaneamente attivi entro la stessa BS, ed owiamente serviti da altrettanti 
beamformer della stessa sezione ricetrasmettitori BTS (Base Station Transceivers). 
Come vedremo in seguito, nei sistemi GSM/DCS la discriminazione tra utile ed 
interferenti cocanale awiene sfruttando una sequenza di training di 26 simboli noti 
che costituisce il midambolo dei burst di trasmissione, posto di aver assegnato 



ad apertura sintetica, in grado cioe di trasmettere e ricevere in direzione dei terminali 



inseguire un generico utente mobile che si muove nel piano azimutale. 
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cocanale. La suddetta sequenza d appositamente prevista pefN^^rap^lfe del 
sincronismo di simbolo e per la stima della risposta all'impulso del canale. Nel caso di 
utenti cocanale della stessa cella, essi sono anche fra loro sincronizzati, 6 quindi 
possibile localizzarli analizzando la sola sequenza di training. 



Adottando le strategic RCS e SCFR, occorre assolutamente evitare che gli 
interferenti cocanale possano peggiorare la quality del segnale demodulato, 
vanificando in tal modo ii vantaggio della maggiore capacita di traffico consentito. Si 
vuole ancora sottolineare che, nei sistemi utilizzanti un approccio del tipo RCS, 
I'interferenza e dovuta agli utenti cocanale delle celle adiacenti, mentre nei sistemi 
che utilizzano I'approccio SCRF I'interferenza e generata sia dagli utenti cocanale 
della stessa cella, che dagli utenti cocanale di celle adiacenti. 

II canale radiomobile e stato a lungo studiato, e diversi modelli sono stati 
proposti compatibilmente con i piu frequenti ambiti trasmissivi, rurale, urbano, etc. Da 
questi studi risulta che, a causa della natura dispersiva spaziale e temporale del 
canale radio, oltre agli interferenti generati dagli utenti cocanale sul percorso diretto, 
esistono anche quelli corrispondenti agli echi di questi ultimi sui percorsi multipli. 

Particolare attenzione deve essere posta ad eventuali echi ritardati del segnale 
utile, poiche in alcuni esempi di arte nota essi sono relegati al ruolo di interferenti 
cocanale, mentre nell'esempio realizzativo secondo la presente invenzione gli stessi 
sono effettivamente riconosciuti come parte integrante del segnale utile. 

Nel seguito concentreremo I'attenzione sulle stazioni radio base che utilizzano 
antenne a schiera per realizzare le tecniche RCS o SCFR, quesfultima anche 



secondo la piu moderna tecnica SDMA, con assegnazione dinamica dei canali, o in 
concomitanza con I'utilizzo di tecniche tradizionali di ripartizione fissa delle frequenze. 

Partiremo anche dall'assunto che e gia stato risolto il problema della 
sincronizzazione, owero noto ii metodo di acquisfzione del sincronismo di portante, e 
di simbolo, e nota la posizione della sequenza di training all'interno del burst. 

ARTE NOTA 

La natura del canale radiomobile e tale da causare una dispersione spaziale 
(piu direzioni di arrivo degli echi) e temporale (ritardo degli echi) dell'energia del 
segnale che giunge a ciascuna delle antenne. In un contesto caratterizzato dalla 
propagazione lungo percorsi multipli, e ragionevole aspettarsi che i diversi percorsi del 
segnale giungano ad una schiera di antenne con DOA different!. Questa caratteristica 
del canale radiomobile pud essere utilizzata per realizzare un'equalizzazione nel 
dominio spaziale, e possibile ad esempio filtrare il solo percorso a massima energia, 
oppure combinare tra loro opportunamente tutti i segnali provenienti dai differenti 
percorsi per cancellare gli interfered cocanale. L'equalizzazione della dispersione 
temporale verra considerata success ivamente. 

Onde poter effettivamente discriminare I'utile dalle diverse tipologie di 
interferenti, occorre dotare la schiera della capacita di modificare dinamicamente la 
propria funzione di direttivita e di inserire degli zeri nel proprio diagramma di 
radiazione nelle direzioni dei segnali interferenti. Nel seguito verranno illustrati alcuni 
esempi realizzativi in tal senso, secondo I'arte nota, evidenziando di volta in volta i 
rispettivi inconvenienti. La descrizione, per semplicita di trattazione e senza nulla 
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togliere alia generality della stessa, 6 limitata ad una array lineare ed uniforme di 
antenne, detto anche ULA (Uniform Linear Array). 

Con riferimento alia fig.1, I'array consiste di M element! uguali allineati ed 
equispaziati. La separazione tra elementi adiacenti e d. Per segnali a banda stretta 
con frequenza della portante fc, si puo definire la corrispondente lunghezza d'onda 
come: 

dove c 6 la velocity di propagazione (c = 3.108 m/s nel vuoto). La spaziatura d fra gli 
elementi puo essere convenientemente espressa in funzione della lunghezza d'onda 
X come: 

d-Ax\ (2) 
A 6 cosi la separazione fra gli elementi in termini di lunghezza d'onda. 

Supponiamo nel seguito che sia presente un'unica sorgente puntiforme, 
sufficientemente distante da rendere assimilabile ad un'onda piana quella incidente 
suirarray. La direzione indicata in fig.1 dairangolo 8 e la direzione di arrivo dell'onda 
piana, detta in seguito DOA (Direction Of Arrival). Si ottiene quindi che il ritardo 
relativo T tra i segnali ricevuti da due elementi adiacenti dell'array 6: 

<*cos(9) 

T = (3) 

C 

II segnale reale riceyuto sul k-esimo element© dell'array 6 allora: 

r k (t)-s(t-T k )-s(t-kz) (4) 
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Si assume che i ritardi relativi ai diversi ricevitori lungo la schiera siano piccoli 
rispetto al reciproco della banda del segnale, cosicche gli stessi si possono 
approssimare con delle opportune variazioni di fase deirinviluppo complesso r k (t) . 



Risulta quindi: 



(5) 



# 



il fattoredi fase2rt/ c x pu6 essere scritto in piu modi: 



2^ x _ 2kf c = 2* ^-cos(e) - 2^008(0) 

C C A- 



(6) 



In definitive se s(t) 6 I'equivalente in banda base del segnale ricevuto su| 
primo elemento dell'array, e definito un vettore a(e) come: 



i(e)- 



i 

e -y23tAcos(e) 

-y23t(A/-l)Acos(6) 



(7) 



possiamo raggruppare i segnali ricevuti sui vari elementi dell'array in un vettore r(r): 

¥(t) ~ a(6)j(r) (8) 

II vettore a(e)e il cosiddetto steering vector, o vettore di risposta dell'array, ad 
un'onda piana incidente con DOA uguale a 6 . 

Nell'arte nota uno dei primi approcci al problema della ricezione e trasmissione 
spazialmente selettive e quello dei cosiddetti sistemi a Switched Beams, owero che 
utilizzano schiere di antenne per formare un numero fisso di fasci orientati in direzioni 
fisse. I fasci cos! ottenuti sono selezionati, sia in trasmissione che in ricezione, in 
modo da utilizzare quel fascio che e meglio orientato nella direzione del terminale 
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mobile Questi sistemi, bench6 di facile realiz^tone, hannjp^ i&i grosso 
inconveniente di non saper discriminare in modo efficace i seg(te|^te^ei^ti 

Con riferimento alia fig.2, notiamo un ricevitore a Switched Beams 
comprendente un array a quattro elementi, connessi agli ingressi di un blocco 
Beamformer, le cui quattro uscite sono rispettivamente connesse agli ingressi di un 
blocco df selezione del fascio ed a quelli di un selettore. L'uscita del blocco di 
selezione del fascio 6 connessa ad un ingresso di controllo del selettore. 

Nel funzionamento il beamformer di fig.2 viene utilizzato per formare un 
diagramma di radiazione composto da un numero fisso B di fasci in direzioni 

prefissate ^ sul piano azimutale. Gli ingressi del beamformer sono rappresentati dal 
vettore f(t) (indicato nel seguito con y(t)) t trascurando per semplicite il segno di 
tilde deirinviluppo complesso. 

II beamformer e costituito da una matrice di coefficient definiti come 




segue: 



w , « e -J2*(~D**»M (9) 



con m = 1, 2, .... M e b = 1, 2, B. II segnale sulla b-esima uscita del beamformer 
sarci dunque pari a: 

y * (0 = ^5(r)e" y2 " (OT " l) ^ cos( * fc) " cos(e) ^ (10) 

m-l 

Se 6 «<j) 6 Tuscita risulta massima. Nel caso ci siano piu segnali, provenienti da 
piu direzioni, considerando la direzione d'arrivo 6 =<|> fr , i contributi dei diversi elementi 
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dell'array si sommano in fase. L'uscita b-esima del beamformer sara dunque 
caratterizzata da un diagram ma di radiazione con un massimo nella direzione ^ . 

La selezione del fascio migliore awiene | secondo un criterio che tende a 
minimizzare la probability d'errore, owero applicando una metodologia MLSE 
(Maximum Likelihoo d Sequence Estimation) di stima della sequenza trasmessa a 
massima verosimiglianza. La suddetta stima e vantaggiosamente compiuta 
applicando il ben noto algoritmo di Viterbi ad un diagramma sequenziale a traliccio, o 
trellis, i cui rami rappresentano tutte le possibili transizioni di stato del ricevitore. 
L'algoritmo calcola, per tutte le possibili sequenze {*„}, la distanza tra il 
corrispondente segnale trasmesso: 

y(<\{x*})=I n x * h ('- nT) 

ed il segnale ricevuto y(t) , e sceglie poi la sequenza {*„} che corrisponde alia 
minima tra queste distanze. Da cio si intuisce che maggiori sono le distanze tra i 
possibili segnali {*„ }) , minore risulta la probability di scegliere durante la 
demodulazione la sequenza sbagliata. Definendo quindi con d^ B il valore minimo 
della distanza fra tutti i possibili y(t\ {x„}), tale parametro risulta essere un indice 

affidabile della probability d'errore. 

A causa della complessita computazionale di questo tipo di approccio, in alcuni 
esempi noti viene anche proposto un altro metodp per la selezione del fascio, basato 
semplicemente sulla qualita della demodulazione della sequenza di training. Fra tutti i 
fasci possibili viene scelto infatti quello che determina il minor numero di errori nella 



clemoclulazione del miclambolo. A parity di errori commessi si utilizza il fascio che, aila 
fine dell'algoritmo di Viterbi applicato al midambolo, da luogo alia distanza minima. 

Un approccio che fornisce all'antenna un maggior livello di intelligenza rispetto 
ai metodi di Switched Beams, consiste essenzialmente nell'esecuzione di un filtraggio 
spaziale sui segnali ricevuti dalle antenne, ad opera di un beamformer che sintetizza 
nel segnale uscente le caratteristiche di direttivita di un opportuno diagram m a di 
radiazione in grado di ridurre I'interferenza cocanale. Questo approccio, nella maggior 
parte degli esempi noti di applicazione, richiede la conoscenza preventiva delle DOA 
di utile ed interferenti cocanale. In letteratura si trovano diversi tipi di algoritmi studiati 
per risolvere il problema della stima delle DOA, avendo a disposizione una schiera di 
antenne riceventi. Ualgoritmo piu rappresentativo 6 quello noto con il nome di MUSIC 
(Multiple Signal Classification), proposto da R. O. Schmidt nell'articolo intitolato 
"Multiple Emitter Location and Signal Parameter Estimation", pubblicato in IEEE 
Trans. Antennas and Propagation, vol. AP-34, No. 3, March 1986, a cui si rimanda 
per una trattazione piu completa. Ualgoritmo necessita dell'ipotesi restrittiva che la 
matrice di covarianza del segnale non sia singolare; in tal caso i segnali sulle antenne 
non sono correlati tra loro e gli steering vector risultano linearmente indipendenti. 

Note le direzioni d'arrivo (e), i = 1, 7 p) bisognera costruire un diagramma di 

radiazione appropriato. Una possibility e quella di definire tanti diagrammi quante 
sono le DOA stimate (p), in base alia tecnica MVDR (Minimum Variance Distortionless 
Response), descritta ad esempio da S. Haykin nel volume intitolato "Adaptive Filter 
Theory", edito da Prentice-Hall, Inc., Englewood Cliffs, N.J., 1991. 
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Tale approccio permette di stimare i coefficient! del beamformer in modo che il 
diagramma di radiazione abbia risposta unitaria in corrispondenza della direzione 
prescelta, minim izzando contemporaneamente Teffetto d'interferenza e del rumore 
proveniente da tutte le altre direzioni. 

Fino a questo punto, quindi, si sono stimate le direzioni di arrive dei vari segnali, 



costruendo per ognuno di essi un apposito diagramma di radiazione. Con la tecnica 
descritta, non c'e modo pero di sapere qual e la direzione di arrivo associata al 
segnale ritenuto utile, e di conseguenza qual e il diagramma di radiazione che e piu 
conveniente utilizzare per effettuare la demodulazione. Si pud quindi pensare di 
selezionare il fascio secondo le modality precedentemente descritte, in modo da 
scegliere quello che permette una migliore demodulazione della sequenza di training. 

Nel caso di stima delle DOA, le argomentazioni di cui sopra sono rigorosamente 
valida per il caso di p segnali incorrelati con p < M, dove M e il numero di elementi 
deirarray, owero per scenari caratterizzati da un basso numero di percorsi multipli, 
come i sistemi di comunicazione mobile in area rurale. 

Nel caso di propagazione in ambienti caratterizzati da numerosi percorsi 
multipli, come ad esempio nelle aree metropolitane, I'estrazione delle DOA pud 
risultare inefficiente, poiche il numero di direzioni di arrivo dei vari segnali e spesso 
superiore alia dimensione deirarray, e anche Pjpotesi di incorrelazione tra i diversi 
segnali sulle antenne 6 difficilmente verificata. Ad ogni modo, almeno nel caso in cui i 
multipercorsi relativi ad uno stesso segnale appartengono tutti ad un settore angolare 
non troppo ampio, le prestazioni delPalgoritmo rimangono spesso soddisfacenti, e 
comunque superiori a quelle garantite dall'approccio a Switched Beams. 
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Per rendere piu efficiente la stima delle direzioni nel caso di segnali ^frelatT, 
pud pensare di utilizzare altri algoritmi di ricerca delle DOA, come il WSF ^ 
Subspace Fitting) proposto da B. Ottersen e M. Viberg neirarticolo intitolato 
of subspace fitting based methods for sensor array processing", a seguito dei Proc. 
ICASSP 89, pp. 2807-2810, Glasgow, Scotland, May 1989. Tale algoritmo, al prezzo 




di una complessita computazionale superiore, riesce a trattare anche sorgenti 
coerenti. I vantaggi vengono comunque fortemente limitati sia, come gik detto, dal 
grande numero di direzioni di arrivo, sia dal fatto che, anche riuscendo a stimare piu 
di una direzione per ogni segnale, bisognerebbe utilizzare un algoritmo per la 
composizione del fascio in grado di sfruttare appieno la conoscenza di tali direzioni. 
Disponendo di uno schema di demodulazione efficiente in grado di gestire i 
multipercorsi, come e quelle basato su un algoritmo MLSE, pud risultare infatti utile 
puntare il fascio verso piu repliche dello stesso segnale, in modo da massimizzare il 
rapporto segnale-disturbo, cosa che pero non e possibile con la gia menzionata 
tecnica MVDR che, come gi& detto, considera come utile una sola direzione. 

Gli algoritmi sopra illustrati che prevedono la stima delle DOA si dimostrano 
inoltre inefficient nel caso di non-stazionarieta durante la trasmissione. Consideriamo 
infatti il caso di un interferente che disturba la trasmissione solo in un semiburst (e 
non durante il midambolo). La direzione d'arrivo di un interferente stimata utilizzando 
un algoritmo di selezione del fascio basato sulla qualita della demodulazione del 
midambolo, potrebbe non considerare I'effetto di tale interferente. In tal caso 
occorrerebbe ricorrere a piu onerosi algoritmi per la selezione, o composizione, del 
fascio basati sulla demodulazione non piu del solo midambolo ma dell'intero burst. 
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Oalio Giurtirv 

Nel sistema radiomobile celluiare, il numero dei segnaii interferenti pud risultare 
a volte maggiore del numero deile antenne che compongono la schiera, in tal caso 
I'ipotesi p < M non 6 piu verificata e I'array pu6 non essere in grado di inserire zeri nel 
diagramma di radiazione in direzione di ogni interferente. Cid significa che i metodi 
basati esclusivamente sulla stima delle DOA perdono ulteriormente di efficacia. 



Tuttavia anche in questo contesto I'array adattativo deve poter essere in grado 
di combinare i segnaii presenti alle antenne per migliorare il rapporto tra il segnale e la 
somma degli interferenti col rumore, rapporto noto in letteratura con I'acronimo di 
SINR, (Signal to Noise and Interferents Ratio). Nei seguenti articoli di J. H. Winters: 

• "Optimum Combining in Digital Mobile Radio with Cochannel Interference", 
pubblicato in IEEE J. Select. Areas Commun., vol. SAC-2, No. 4, July 1984; 

• 'Optimum Combining for Indoor Radio Systems with Multiple Users", pubblicato in 
IEEE Trans. Commun., vol. COM-35, No. 11, November 1987, viene sviluppato il 
concetto di combinatore ottimo. 

E questa una tecnica alternativa alia stima delle DOA, che si basa sulla 

conoscenza "a priori" da parte del ricevitore di una sequenza di simboli noti 
(sequenza di training), presente nei singoli burst di trasmissione per scopi di 
sincronizzazione. La sequenza viene ulteriormente sfruttata nel processo di 

filtraggio spaziale mediante il quale vengono opportunamente combinati i segnaii. 
ricevuti dalla schiera di antenne per massimizzare il SINR. 



Tale approccio consente di sfruttare la diversity spaziale del segnale, oltre che 
per ridurre la potenza relativa degli interferenti cocanale, anche per affrontare il fading 
alia Rayleigh del segnale utile. 

In fig. 3 viene mostrato lo schema a blocchi di un combinatore ottimo a M 
element!. Come si pud notare, il combinatore ha M ingressi per altrettanti segnali di 



banda base y 1 (n) 9 ... 9 y.(n\...,y M (n) 9 ricevuti da corrispondenti M elementi di un 
array (non visibile in figura). I suddetti segnali vengono moltiplicati per dei rispettivi 
pesi w 1 ,...,w i ,...,w w , in genere complessi. I segnali pesati pervengono agli ingressi 
di un primo sommatore, alia cui uscita e presente un segnale y oui (n) che fuoriesce 
dal combinatore e perviene inoltre, con segno negato, ad un ingresso di un secondo 
sommatore, al cui altro ingresso e presente un segnale x n . All'uscita del secondo 
sommatore 6 presente un segnale errore e n che perviene ad un blocco WEIGHT 
ADJUSTMENT avente altri ingressi per i segnali y 1 (n\... 9 y i (n) 9 ...,y M (n) . Alle uscite 
di detto blocco sono present! i pesi w 19 ... 9 w. 9 ... 9 w M - 

Nel funzionamento, il segnale y out (n) pud essere espresso come: 

Indicando con w-[yv p »v 2 , w tf J il vettore dei pesi, e con 

y(n)^[y l (n\y 2 (n\.^.,y M (n)J il vettore dei segnali ricevuti, ed adottando una 
notazione matriciale, possiamo definire la somma pesata dei segnali provenienti dalle 
antenne come: 
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^W-w T y(«) (13) 
La differenza tra il segnale desiderate x nPI corrispondente ai campioni della 
sequenza di training, ed il segnale y out {n) in uscita dal combinatore definisce il 
segnale errore: 

. "Mm) (") 



Questo segnale viene usato come segnale di controllo nel processo di 
adattamento dei pesi del combinatore. 
© Lo scopo del processo di adattamento, § quello di trovare un vettore di pesi w, 

tale da fornire in uscita, ad ogni istante, un segnale uguale a quello desiderato x n . 
Quando i segnali in ingresso possono essere visti come variabili stocastiche 
stazionarie, il blocco WEIGHT ADJUSTMENT cerca il vettore dei pesi che minimizza 
una funzione costo corrispondente all'errore quadratico medio MSE (Mean Square 
Error), definito come segue: 

e 2 (w) = E^e m \ 2 ] (15) 

dove £[■••] indica il valore medio dell'argomento tra parerrtesi. 

Sostituendo la (14) nella (15) possiamo ottenere I'insieme dei pesi che 
minimizza I'MSE: 



e 2 (w) = E 
= E 



\ x n\ 2 +w'y(n)y'Ww-2w\y(n) 

21 - ' - (16) 
kl J+ WR w w-2w'R^ 



dove la matrice di covarianza spaziale R w dei segnali in ingresso alle antenne pud 
essere espressa come: 
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R yy * E[yb*y(n)] 



e la matrice di cross-covarianza R tra gli M segnali in ingress 



segnale desiderata, pu6 essere espressa come: 




ne ed il 



(18) 



L'errore quadratico medio definito nella (16) e una funzione quadratica dei pesi 
w. Le componenti del gradiente della funzione MSE, sono le derivate parziali dell'MSE 
rispetto ai valori dei pesi. Derivando la (16) rispetto a w, otteniamo il gradiente come 
funzione lineare dei pesi: 



v4„r] - 



2R„w-2R^ 



(19) 



Quando la scelta dei pesi e ottimizzata, il gradiente risulta nullo, quindi: 



R 



(20) 



II vettore dei pesi w e dunque quello chefornisce il minimo errore quadratico 

medio. Uequazione (20) b detta equazione di Wiener-Hopf. Per un'analisi piu 
dettagliata dei metodi di ricerca del vettore dei pesi ottimo, si rimanda airarticolo 
intitolato "Adaptive Antenna Systems autori: B. Widrow, P. E. Mantey, L J. Griffiths 
e BB. B. Goode, pubblicato negli IEEE Proceedings, vol. 55, No. 12, December 1967. 
Considerando ora I'equalizzazione temporale, 6 opportune introdurre le 

seguenti considerazioni: sia {x n } una sequenza di N simboli trasmessi su un canale 
con lunghezza di dispersione L, owero la cui risposta impulsiva sia significativa entro 
un tempo LT, essendo T il tempo di simbolo. Consideriamo una schiera di antenne 
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composta da M element!. II vettore dei segnali ricevuti sulla schiera di antenne 



y(r) -[*('), J>2('), 



\,.... 9 y M (t)T 6 rappresentato dalla convoluzione di un vettore h(t) 



con la sequenza trasmessa 



y(0- 2*Ji(f-»r) + n(0 



(21) 



dove gli elementi del vettore h(t -nT) sono i campioni della risposta impulsiva del 
canale sulle varie antenne. II vettore h(/) puo essere scritto come: 



essendo: 

• N e ii numero dei cammini multipli, 

• A l Tattenuazione di tratta deiri-esimo cammino, 

• G ( . la direzione d'arrivo, 

• h 0 la convoluzione dell'impulso di modulazione e della risposta impulsiva del filtro 
di ricezione, 

• a(e ; ) 6 lo steering vector, o vettore della risposta dell'array in direzione 9 ; , definito 



II vettore n(r) rappresenta il rumore additiyo, che considereremo 
temporalmente e spaziaimente bianco, e che soddisfa quindi alle seguenti condizioni: 



• 




N 



(22) 



nella (7). 




(23) 
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In gran parte degli esempi noti di realizzazione di un equalizzatore ottimo per la 
compensazione delle distorsioni introdotte da! canale sulle caratteristiche temporal! 
della sequenza trasmessa, I'equalizzatore e costituito da un filtro adattato alia risposta 
impulsiva del canale, seguitp da un decisore a massima verosimiglianza (MLSE) 
basato sull'algoritmo di Viterbi. Sia per la realizzazione del filtro adattato, che per 



quella del decisore MLSE, 6 quindi indispensabile conoscere la risposta impulsiva del 
canale, in genere ignota. Sfruttando la conoscenza del midambolo, I'equalizzatore 
stima la risposta ail'impulso del canale, per poi calcolare i coefficient! del filtro adattato 
alia risposta impulsiva del canale e poter quindi stimare la sequenza trasmessa, 
vantaggiosamente utilizzando ralgoritmo di Viterbi. La possibility di effettuare una 
stima del vettore h(r) rende quindi superfluo il calcolo esplicito dello stesso secondo 
la (22), che presuppone la stima delle direzioni d'arrivo. 

Gli algoritmi finora menzionati, affrontano in due fasi di equalizzazione disgiunte 
il problema della stima della sequenza trasmessa, e precisamente: 

• equalizzazione spaziale per Pannullamento degli interferenti; 

• equalizzazione temporale del segnale filtrate spazialmente. 

Un'equalizzazione disgiunta consente pero di combinare spazialmente solo gli 
echi del segnale che non risultino significativamente ritardati. Gli echi del segnale 
ritardati sono infatti considerati come interferenti cocanale. Tale approccio pu6 pero 
risultare in una riduzione del rapporto tra segnale e rumore all'ingresso 
deirequalizzatore temporale, come verr& visto confrontando tra loro le figg.4 e 7, ed in 
una minore robustezza dell'algoritmo di beamforming a fronte di particolari 
distribuzioni azimutali degli echi, come verr& visto confrontando tra loro le figg.5 e 8. 



Italtel sp. 

Alia luce di quanto appena detto, prendiamo in considerazione la funzione 
costo (15) del combinatore di Wiener-Hopf, discussa in precedenza: 




w.y.(n)-.x H 



2 1 



(15) 



dove con w i si indica il coefficiente con cui si pesa il segnale ricevuto alia i-esima 



antenna. ;Questo algoritmo massimizza il rapporto tra il segnale ed il rum ore, dove per 
segnale si intendono tutti quegli echi che risultano correlati con la sequenza di training 

{*„}. Se il canale introduce dei ritardi superiori al tempo di simbolo, Is 3J\is per il 
GSM, gli echi cosl ritardati risultano decorrelati con la sequenza di training e vengono 
trattati dal combinatore come interferenti cocanale. 

In fig.4 e rappresentato il diagramma di radiazione sul piano azimutale di una 
schiera di antenne, in cui i coefficient! di beamforming sono caicolati con I'algoritmo di 
Wiener-Hopf applicato al funzionale (15), in presenza di un segnale utile generato da 
un mobile MSd, con relativo eco ritardato di 5 \is t e di un interferente cocanale 
generato da un mobile MSi. Come si pud notare, il diagramma di radiazione presenta 
un massimo nella direzione dell'utile (90°) , un primo zero nella direzione 
dell'interferente cocanale (120°), ed un secondo zero nella direzione dell'eco del 
segnale utile (60°). A causa del secondo zero si perde la possibility di sfruttare anche 
I'energia dell'eco dei segnale utile nel processo di stima della sequenza trasmessa. 

In fig. 5 6 rappresentato il diagramma di radiazione sul piano azimutale della 
schiera di antenne utilizzata in uno scenario che differisce dal precedente solo per la 
presenza di una replica dell'interferente cocanale che viene a trovarsi sulia stessa 



direzione di arrivo (90°) del segnale principale. ComeW^pJa^ notar^yil^dfagramma di 




radiazione presenta ancora i due zeri nelle direzioni pf&oedl&jj&r nspettivamente, 
delllnterferente (120°) e dell'eco del segnale utile (60°) , ma in questo nuovo scenario 
non 6 piu presente il massimo nella direzione dell'utile (90°). Questo com portamento 
si spiega, considerando un sistema di sincronismo tradizionale, in quanto il 




combinatbre vede un unico segnale utile, distinto dall'interferente, nella comune 
direzione (90°), e tenta comunque di ridurre i'effetto delllnterferente in quella 
direzione riducendo il guadagno. Cid com porta un'ulteriore riduzione dell'energia del 
segnale utile presente nel processo di stima della sequenza trasmessa. II combinatore 
perde quindi ['opportunity di trasferire il lobo ,di massimo guadagno nella sola 
direzione dell'eco del segnale utile (60°) , riuscendo in tal modo ad annullare I'effetto 
dell'eco delllnterferente avente la stessa direzione dell'utile. 

In ambito radiomobile sono noti degli algoritmi piu recenti di quelli finora illustrati 
che cercano di owiare agli inconvenient! evidenziati dalla discussione di figg.4 e 5. I 
suddetti algoritmi sfruttano ancora il concetto di minim izzazione dell'errore quadratico 
medio espresso da un opportuno funzionale, come aweniva per il combinatore 
ottimo, ma, a differenza di questi, I 'inform azione sul canale che entra nell'espressione 
del suddetto funzionale e ora fornita da un filtro cosiddetto spazio-temporale per il 
fatto di doversi adattare alia risposta complessiva del canale. 

Un esempio in tal senso 6 costituito dalla domanda di brevetto internazionale 
W0 97 08849 A, pubblicata il 6/3/97 a nome THOMSON CSF, intitolata: "METODO E 
DISPOSITIVO PER IL MULTIPLEXING-DEMULTIPLEXING SPAZIALE DI SEGNALI 



m 



RADIOELETTRICI PER SISTEMI RADIOMOBILE SDMA", con priority francese del 
22-8-95. 

Nella suddetta domanda viene descritta una stazione radio base che utilizza in 
maniera "intelligente" una schiera di antenne a cui e accoppiato un ricetrasmettitore 
multivia dotato di mezzi di beamforming particolarmente idonei a rendere possibile 



I'assegnazione di un medesimo canale frequenziale-temporale a piu mobili nella 
stessa cella, sempre che i relativi segnali abbiano direzioni d'arrivo sufficientemente 
decorrelate L'invenzione, che sostanzialmente rivendica un'applicazione della tecnica 
SDMA in ambito radiomobile, presenta due variant! realizzative che si distinguono per 
la diversa tipologia dei mezzi di beamforming utilizzati in ricezione. Piu precisamente, 
la prima variante e riferita a mezzi che necessitano di un blocco funzionale per la 
stima degli steering vector, che vengono quindi passati ad un filtro spaziale di 
beamforming vero e proprio. I mezzi di beamforming della prima variante sono 
riconducibili a quelli g\h menzionati per la stima delle DOA e risultano molto meno 
pertinenti nei confront! della presente invenzione rispetto ai mezzi di beamforming 
previsti nella seconda variante, pertanto la prima variante non verr& esaminata. Giova 
tuttavia ricordare che le tecniche utilizzate nella prima variante sono dette "trasparenti" 
perche non necessitano di alcuna conoscenza^a priori sulla sequenza di simboli 
ricevuta. 

I mezzi di beamforming previsti nella seconda variante realizzativa della 
citazione THOMSON utilizzano, invece, ia conoscenza a priori della sequenza di 
training, come appunto awiene nell'invenzione che sar& descritta. Questi secondi 
mezzi comprendono un filtro spazio-temporale. L'ottimizzazione della risposta 



complessiva del canale viene compiuta mediante due algoritmi tra loro alternative 
rispettivamente denominati FAST-BBT-DFE e FAST-DFE, dove: 

• FAST 6 un acronimo di Filtro Adattato Spazio-Temporale, 

• BBT e un acronimo di Bruit Blanc Temporellement, con terminologia francese, e 

• DFE di Digital Feedback Equalizer, con terminologia inglese. 

Riportiamo di seguito quanto e testualmente detto in merito ai vantaggi di questi 
algoritmi: 

• "II FAST-BBT-DFE permette di rimettere in fase I'energia dell'insieme delle 
traiettorie del segnale utile e di reiettare in modo sub-ottimale i segnali interferenti 
da percorsi multipli. Questo modo di procedere comporta un guadagno in termini di 
rapporto segnale rum ore sul simbolo corrente in confronto ad un algoritmo di 
filtraggio solo spaziale (FAS-DFE), infatti I'algoritmo presenta le stesse capacita di 
reiezione dei segnali interferenti con la differenza che gli echi dell'utile diversi dalla 
traiettoria principale non sono piu assimilati a dei segnali interferenti, cosicche si 
liberano altrettanti gradi di liberta per prendere in considerazione altri interferenti. 

• II FAST-DFE permette di trattare in modo quasi-ottimale un segnale con interferenti 
proveniente da un array, in quanto esso realizza non solo la messa in fase degli 
echi dell'utile ma anche quella degli echi degli interferenti prima della loro 
reiezione. Questo modo di procedere consente di eliminare una sorgente 
interferente utilizzando un solo grado di liberty spaziale qualunque sia il numero di 
echi secondo i quali questa sorgente 6 ricevuta dairarray". 

Occorre tuttavia precisare che le suddette operazioni di "messa in fase" o 
"riallineamento" sono soltanto il risultato degli algoritmi preposti e non costituiscono 
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quindi passi autonomi del processo di ottimizzazione del canale, poich6 altrimenti 
richiederebbero una stima separata delle DOA. Inoltre, il blocco di stima delle DOA 
resta anche nella seconda variante per gli scopi tipici della tecnica SDMA, owero per 
poter compilare una tabella di correlazione intermobili delle direzioni di arrivo dei 
segnali alle antenne, da utilizzarsi nel multiplexing spaziale per il beamforming di 



trasmissione. 

L'algoritmo FAST- B BT- D F E , illustrato con riferimento alia fig. 11 della citazione, 
parte dalla formulazione teorica di un filtro ottimo (valendo le ipotesi restrittive del 
caso) per ciascuna delle N vie facenti capo alle N antenne. II suddetto filtro non 
agisce direttamente sul segnale ricevuto, ma bensi sullo stesso prima sottoposto ad 
uno sbiancamento spaziale della relativa componente di rumore+interferente. Per 
compiere il filtraggio ottimo 6 necessario conoscere la matrice R b di autocorrelazione 
di rumore+interferente sulle N antenne e la matrice G di risposta compiessiva del 
canale. Vengono alio scopo introdotte delle ipotesi restrittive sulla natura della matrice 
R che consentono di stimare separatamente le due matrici e completare quindi il 
filtraggio proposto all'inizio del procedimento. Dopodiche il segnale filtrato viene fatto 
pervenire ad un organo di decisione di tipo DFE. 

La principale ipotesi restrittiva consiste nel considerare la matrice R t di natura 
impulsiva, il che comporta degli interferenti tra loro completamente incorrelati, e quindi 
equiparabili ad un rumore bianco nella banda del segnale utile. In altri termini, gli 
interferenti sono considerati privi di echi, o meglio, ciascun eco e considerate come 



• 



esso stesso un autonomo interferente. Sotto queste ipotesi la convergenza 
deiralgoritmo 6 detto essere abbastanza rapida. 

E del tutto evidente che I'applicazione di un tale algoritmo non 6 affatto 
consigliabile in ambito radiomobile dove, specie nelle celle urbane di 
dimensioni, lo scenario tipico 6 al contrario caratterizzato da interferenti jb©c&nale 



corredati da un ventaglio di relativi ecni da percorsi multipll. thfatto di non co 
questa realty com porta lo svantaggio di non poter ottimizzare in modo effettivam 
credibile il rapporto segnale-rumore complessivo, : oppure, qualora venga utilizzata la 
tecnica SCFR (in modality SDMA) di non poter separare in modo ottimale sul piano 
azimutale i diversi utenti mobili cocanale. 

Ualgoritmo FAST-DFE, illustrato con riferimento alia fig. 12 della citazione, non 
pone invece alcuha limitazione sulla natura del rumore, comprensivo questi degli 
interferenti cocanale, trattando con cio il canale in maniera quasi ottimale. Pure in 
questo caso la struttura e costituita da un filtro spazio-temporale associato ad un 
organo di decisione di tipo DFE, anche se e detto possibile I'utilizzo del Viterbi, senza 
perd illustrare le modality di un tale utilizzo. II filtro spazio-temporale 6 costituito da N 
filtri temporali W AI lunghi N HR , uno per ogni via n-esima che si diparte dalle N 
antenne. ^ inoltre presente un ulteriore filtro temporale HR, pure lungo N m posto 
sulla via di retroazione deirequalizzatore DFE. Come § noto dalla teoria degli 
equalizzatori, un DFE, prima di sottoporre i simboli ricevuti ad un decisore a soglia, 
cerca di rimuovere HSI dalla sequenza dei relativi campioni, predisponendo alio scopo 
un filtro cancellatore che agisce in retroazione sui campioni postcursori della risposta 
impulsiva, e di un equalizzatore prelim inare che agisce in avanti (forward) sui 




# 



26 




campioni precursori. Nello schema di fig.12 il fiitro HR costituisce appunto il 
cancellatore, mentre I'insieme dei filtri W*, i cui segnali filtrati vengono sommati tra 
loro a costituire un unico segnale somma, sono equiparabili all'equalizzatore 
preliminare. II segnale ripulito dall'ISI, indicato nel testo con i(n), e ottenuto dalla 
differenza tra il suddetto segnale somma ed il segnale filtrato da HR. II segnale z(n) 



attraversa il decisore a soglia che fornisce una stima z(n) della sequenza di dati 



trasmessa. 

L'aggiornamento dei pesi dell'intera struttura filtrante viene effettuato 
minimizzando I'errore quadratico medio MSE espresso dal seguente funzionale: 

£ =E[|z(t)-z(tf] (24) 
II funzionale (24) nel caso di canale stazionario viene stimato come segue: 

i =l2|z(i)-z(if (25) 



mentre per canale non stazionario vale la seguente espressione: 

£ = Itx ( "- i) |z(i)-z(i)| 2 =X^(/»-l) + |z(i)-z(if (26) 

dove X k un fattore di oblio scelto tra zero e uno in funzione della non stazionarieta 
del canale. 

In presenza della sequenza di training, questa sostituisce z(i) nelle (25) e (26). 

Manca nella descrizione della seconda variante una chiara indicazione di come 
agire sui diversi filtri per minimizzare I'MSE, a scopo di descrizione sono stati citati gli 
acronimi dei principali algoritmi noti di minimizzazione ma non viene fornito un valido 
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esempio di applicazione al caso specifico spazio-temporale. Manca inoltre" *° 
I'indicazione di un vincolo da imporre sull'energia del segnale ricevuto, al fine di 
evitare nella minim izzazione della (24) la possibility di convergenza verso una classe 
di soluzioni banali. Posto che nella citazione THOMSON gli esempi realizzativi con 
filtraggio spazio-temporale sono stati promossi al rango di variante dotata di 



autonoma originalita, e che I'applicazione aegli algoritmi anzidetti non sembra nota al 
caso spazio-temporale, la richiedente dubita che la semplice citazione degli stessi 
possa mettere in grado il tecnico del ramo di poter riprodurre I'invenzione rivendicata 
nella variante senza esercitare autonoma attivit& inventiva, 

ANAUSI DEL PROBLEMA TECNICO 

Quelle che appare fin da subito evidente 6 una certa qual complessit& 
dell'intera struttura di fig. 12 della citazione THOMSON. Analizzando infatti 
I'espressione del funzionale da minim izzare, nonch6 la stessa figura, e come h detto 
d'altronde chiaramente nel testo, si vede che I'ottimizzazione awiene congiuntamente 
sul filtro spazio-temporale W e sul filtro HR deirequalizzatore DFE. Occorre quindi 
aggiornare un numero di (N + l)xN HR coefficient! indipendenti ad ogni intervallo di 
simbolo. Uelevato numero di coefficient! da aggiornare com porta una certa qual 
lentezza nella convergenza della stima e quindi la necessity di sequenze di training 
lunghe, fatto questo che viene autonomamente evidenziato in una tabella dei 
vantaggi ed inconvenienti dell'algoritmo FAST-DFE. 

Volendo utilizzare il suddetto equalizzatore spazio-temporale in ambito GSM, 

dove le sequenze di training previste dallo standard hanno categoricamente una 

» 

lunghezza di 26 simboli, e possono quindi essere considerate brevi, sorge il problema 
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della convergenza della stima entro la finestra temporale di una tale sequenza. In tal 
caso, infatti, la stima deli'errore £ subirebbe un ulteriore rallentamento dovuto alia 
mancanza di precisione nella valutazione dell'errore stesso, dovendo alio scopo il 
DFE utilizzare delle stime z(n) in sostituzione dei campioni non piu disponibili della 
sequenza di training. Esiste quindi la concreta possibility che una parte consistente di 
simboli dell'intero burst di trasmissione venga erroneamente stimata, causando con 
cid un degrado della qualita del segnale demodulato. Nel caso poi di sistemi 
trasmissivi che adottino dei burst corti, tra cui potrebbe rientrare lo stesso GSM con un 
burst di solo 148 simboli, si profila addirittura if pericolo di una mancata convergenza 
entro I'intera durata dello stesso burst. 

L'inconveniente messo in luce e gia di per se sufficiente a sconsigliare I'utilizzo 
di un siffatto algoritmo in ambito GSM, ma ulteriori considerazioni in merito verranno 
fatte per evidenziare dei problemi tecnici conseguenti e di non poco conto. 

Un secondo problem a che pud nascere dalla lentezza nella convergenza della 
stima e quello della scarsa capacita dell'equalizzatore di fronteggiare le non 
stazionariete del canale, dovute principalmente ad interferenti che si attivano al di 
fuori della finestra temporale della sequenza di training, e/o al caso in cui il canale 
non si possa piu considerare tempo-invariante, come ad esempio awiene quando e 
sensibile I'effetto doppler dovuto all'alta velocita dei mobili. Gli inconvenienti sopra 
evidenziati diventano ancora piu rimarchevoli in presenza di interferenti 
non-stazionari improwisi e di forte intensita (c/J<-10dB). E questa una 
condizione che pud verificarsi frequentemente in ambito radiomobile urbano, ad 
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esempio, quando per via di una cosiddetta "corner attenuation" la potenza degli 
interferenti viene improwisamente a sovrastare quella dell'utile. Specie in quesfultimo 
contesto occorrono soluzioni ben piu potenti di quella indicata nella (26), consistente 
nel pesare ciascuna iterazione della stima con il fattore di oblio X ^ 1. 

Un ultimo inconveniente, originato dalle stesse cause dei precedent!, consiste 




nella diffieolte di realizzare materialmente un ricevitore in grado di svolgere in un solo 
tempo di simbolo tutte le operazioni aritmetiche necessarie airadattamento 
gran numero di coefficienti della struttura filtrante spazio-temporale. 

SOMMARIO DELLI N VENZI ON E 
La presente invenzione si prefigge di superare gli inconvenienti suddel 
come scopo principale quello di indicare un procedimento di equalizzazione spaziale e 
temporale che sia particolarmente rapido nel convergere verso delle prestazioni 
ottimali, e che pud quindi essere utilmente impiegato in stazioni radio base di un 
sistema radiomobile cellulare, ad esempio GSM, dotate di antenne a schiera 

cosiddette intelligent!. 

Secondo scopo, concomitante al primo, e quello di indicare le modality di come 
estendere il procedimento di equalizzazione al di fuori del midambolo, in presenza di 
un canale non stazionario, specie se affetto da interferenti improwisi e di forte 
intensity, garantendo la piu breve durata possibile del transitorio entro cui 
inevitabilmente si manifesta una perdita di qualita del segnale demodulato, e la 
limitazione di tale perdita a valori contenuti. 
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Terzo scopo d quello di indicare un metodo numerico iterativo per 
raggiornamento di una stima dei parametri di equalizzazione ottimi, che sia nello 
stesso tempo preciso e di rapida esecuzione. 

Quarto scopo e quello di indicare le modalita realizzative mediante array sistolici 
di un equalizzatore basato sul metodo numerico iterativo appena menzionato. 

Per conseguire tali scopi la presente invenzione ha per oggetto un metodo di 
equalizzazione per la cancellazione di interferenti isofrequenziali in segnali ricevuti da 
antenne a schiera di stazioni radio base di un sistema radiomobile cellulare. Detti 
segnali essendo ricevuti sotto forma di burst di simboli di un burst di trasmissione, 
comprendente una sequenza di training nota al ricevitore, detti burst ricevuti essendo 
demodulati, digitalizzati, memorizzati, e sincronizzati mediante discriminazione della 

sequenza di training. 

II principale problema tecnico evidenziato nell'arte nota viene risolto dal metodo 
di equalizzazione della presente invenzione, mediante le seguenti fasi operative: 

a) filtraggio spaziale dei segnali provenienti dalle antenne, utilizzando un filtro 
spaziale che combina tra loro sulla base di opportuni pesi dei simboli 
corrispondenti nei diversi burst memorizzati; 

b) filtraggio temporale del segnale generate da detto filtro spaziale, utilizzando un 
filtro adattato a campioni della risposta impulsiva del canale trasmissivo; 

c) ottimizzazione congiunta di detti pesi e campioni sulla base dei simboli di detta 
sequenza di training e dei detti segnali provenienti dalle antenne present! nei burst 
memorizzati; ' 



d) stima dei simboli del burst di trasmissione corrispondenti alia sequenza uscente dal 
filtro temporale adattato; 

e) aggiornamento dei valori ottimizzati di detti pesi e campioni, al di fuori della finestra 
temporale della sequenza di training, mediante ripetizione delle precedent! fasi 
a)...d), sostituendo nella fase c) i simboli della sequenza di training con dei simboli 

= stimati* nella fase d), rendendo con cid detta equalizzazione adattativa, come 

anche descritto nella prima rivendicazione. 

Vantaggiosamente I'organo decisore 6 un processore di Viterbi che sfrutta i 
detti campioni ottimizzati della risposta impulsiva del canale per calcolare, in modo 
noto, delle metriche di ramo di un trellis, utilizzate per stimare a massima 
verosimiglianza la sequenza di simboli trasmessa. Uorgano decisore potrebbe anche 
essere un DFE, ma questo comporterebbe una perdita di prestazioni. 

Detta ottimizzazione congiunta dei parametri di equalizzazione & conseguita 
minimizzando un funzionale corrispondente all'errore quadratico medio che nasce dal 
confronto tra corrispondenti campioni /x-esimi di due sequenze, di cui, una prima e 
ottenuta sottoponendo i simboli della sequenza di training ad un filtro temporale che 
stima i detti campioni della risposta impulsiva del canale trasmissivo, ed una seconda 
sequenza 6 ottenuta sottoponendo i campioni rc-esimi dei detti burst memorizzati ad 

un filtro spaziale che stima i detti pesi. 

Vantaggiosamente, la suddetta stima congiunta dei parametri ottimali di 
equalizzazione viene effettuata in modo iterativo su grandezze numeriche complesse 
espresse in forma algebrica matriciale e rappresentative degli inviluppi di rispettivi 
segnali in banda base. 



II metodo di equalizzazione oggetto della presente invenzione possiede, rispetto 
all'arte nota piu pertinente costituita dal metodo FAST-DFE, il vantaggio di una piu 
rapida convergenza, da cui deriva la capacita di meglio fronteggiare scenari 
trasmissivi particolarmente critici. Questo vantaggio si dimostra fondamentale in 
ambito GSM ed in sistemi che adottino burst di trasmissione particolarmente corti. 

I vantaggi evidenziati derivano essenzialmente dal ratto di eseguire in due tasi 
successive il filtraggio complessivo, una prima fase di filtraggio solo spaziale, ed una 
seconda fase solo temporale,' ma di ottimizzare congiuntamente i coefficient! di 
entrambi i filtri (per quanto riguarda i coefficient! del filtro temporale adattato, gli stessi 
sono indirettamente ottimizzati attraverso I'ottimizzazione dei campioni della risposta 
impulsiva). Viene con cid a determinarsi una reciproca interazione tra i due tipi di 
filtraggio, che spinge verso I'ottimo per effetto del vincolo imposto dal criterio MMSE 
(Minimun-MSE) in concomitanza con i dati certi della sequenza di training. II risultato 
di una tale interazione 6 quello di favorire il passaggio attraverso la cascata dei due 
filtri della massima energia possibile proveniente dal segnale utile, comprensivo quindi 
dei suoi echi, e di reiettare il maggior numero possibile di interferenti cocanale, a cui 
devono quindi essere ricondotti i relativi echi. Quello che in effetti viene compiuto e 
una risincronizzazione automatica degli echi, valida sia per Tutile che per gli 
interferenti, che opera in sinergismo con la discriminazione della direzione dei diversi 
segnali sul piano azimutale. Le caratteristiche funzionali del procedimento in oggetto 
Ip distinguono in maniera originate rispetto ai procedimenti dell'arte nota piO 
pertinente, dove, come si e visto, b prevista un'unica fase di filtraggio 
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spazio-temporale per modellare completamente il canale, a cui pud seguire 
un'ulteriore fase di filtraggio soio temporale nel caso di utilizzo di un DFE. 

Grazie alle sue caratteristiche funzionali, il procedimento secondo ('invenzione 
manipola un numero di coefficient considerevoimente limitato, pari alia somma M+N 
dei coefficienti di due filtri FIR, questo favorisce senza dubbio la velocita di 




o 

per 



convergepza della stima. 

Ulteriori vantaggi relativi ai metodi di calcolo ed al canale non ^Szibnario 

fa 



verranno di volta in volta elencati procedendo nella descrjzione. 

Occorre inoltre far notare che, malgrado la similitudine tra la form 
funzionale (26) e quello utilizzato nella presente invenzione, che sara in seguito 
illustrato e di cui b stata comunque gia data una chiara rappresentazione, c'd in realta 
una diversity nei due termini della differenza che costituisce Pargomento del modulo. 
Nel caso deH'arte nota f infatti, Tazione correttiva del secondo termine solo temporale 
risulta meno efficace nel modificare nel senso deirottimo la componente spaziale, 
rispetto all'azione del termine corrispondente nel funzionale della presente 
invenzione, perche una componente solo spaziale in realta non esiste, essendo il 
primo termine di natura spazio-temporale, ed 6 questa una causa ulteriore che 
rallenta la convergenza 

Da ultimo, il termine "stima congiunta" assume nella presente domanda un 
significato effettivo supportato dalla descrizione, distinguendosi con cid da un 
significato sottinteso che la stessa dizione ha neiriilustrazione deiralgoritmo 
FAST-DFE. La modality esecutiva di una siffatta stima sar& fornita in relazione alia 
descrizione di un metodo innovativo di minimizzazione del funzionale errore, che sark 
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oggetto di relative rivendicazioni dipendenti. Giova anticipare che detto metodo e 
fondato sull'aggiornamento iterativo di una matrice strutturata ad opera di un vettore 
aggiornamento che include i segnali dei burst memorizzati ed i simboli della sequenza 
di training, e dell'utilizzo della matrice aggiornata per il calcolo di detti pesi e campioni 
ottimizzati al tempo di simbolo o-esimo. 

Nel. metodo oggetto principale d'invenzione la stima congiuma dei paramein 
ottimali di equaJizzazione si esaurisce all'ultimo campione della sequenza di training, 
generalmente inclusa in un midambolo del burst di trasmissione. I 26 campioni del 
midambolo GSM (dei 148 disponibili) sono sufficienti a garantire la convergenza della 
stima anche in presenza di segnali ricevuti di mediocre qualita. I parametri ottimali di 
equalizzazione restano poi disponibili per la demodulazione della parte restante del 
burst memorizzato, ed il segnale demodulato non perde di qualita, purche i canali 
sianotempo-invarianti e gli interferenti abbiano caratteristiche stazionarie. 

II metodo dell'invenzione risolve nella fase adattativa e) il problema di una 
corretta equalizzazione pur in presenza di canali con caratteristiche tempo-varianti e/o 
interferenti non-stazionari, aggiornando la suddetta stima dei parametri ottimali di 
equalizzazione ad ogni nuovo simbolo ricevuto. Operativamente, occorre alio scopo 
poter disporre sia di un vettore dei segnali d'ingresso alle antenne, che di un vettore 
dei dati trasmessi. Mentre il primo vettore e sempre disponibile, il secondo vettore e 
completamente noto solo durante il midambolo, ove corrisponde alia sequenza di 
training memorizzata nel ricevitore. Si vede quindi la necessita, una volta esaurito il 
midambolo, di sostituire la sequenza di training con i simboli della sequenza stimata, 
tenendo presente che i simboli gia decisi ad un istante n, altri non sono che quelli 
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associati ad un percorso selezionato nei trellis fino all'istante n-esimo. In 
concomitanza viene inviata al processore di Viterbi una sequenza aggiornata di 
cam pioni, ottimali della risposta impulsiva del canale, per I'aggiornamento dei valori 
delle metriche di ramo. Con cio i vantaggi prima menzionati sono pienamente 
apprezzabili lungo tutto il burst di trasmissione. 




Questo modo di procedere introduce inevitabilmente, a causa del tempo di 
latenza deH'algoritmo di Viterbi che si somma al ritardo dovuto al filtro temporale 
adattato, un ritardo eccessivo nella stima dei parametri ottimali di equalizzazione, 
impedendo con cio di inseguire variazioni veloci delle caratteristiche del canale. 

Uinconveniente viene evitato forzando I'algoritmp di Viterbi a decidere gi& dopo 
pochi passi (tipicamente 4, o 5) qual e la sequenza di dati piu probabile, che viene 
presa al tempo di simbolo attuale risalendo nel trellis il percorso a massima 
verosimiglianza, mentre le decision! definitive per la decodifica MLSE vengono invece 
prese alia fine del burst, in modo che siano ii piu possibile affidabili. 

Poiche le decisioni a breve termine possono risultare errate, si pub generare un 
fenomeno di propagazione degli errori, che pu6 durare fintanto che i dati erronei 
vengono riutilizzati per la stima dei parametri ottimali di equalizzazione, owero per un 
numero di tempi di simbolo pari al ritardo intrinseco con cui si rendono disponibili le 
decisioni a breve termine. 

II metodo di equalizzazione oggetto d'invenzione elimina anche questo secondo 
inconveniente ricopiando in blocco, ad ogni tempo di simbolo, gli L simboli necessari 
airaggiornamento dei parametri di equalizzazione. E del tutto evidente che 
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raggiornamento dei parametri ottimali di equalizzazione altri non e che una nuova 
stima congiunta degli stessi. 

Vantaggiosamente, e possibile procedere al suddetto aggiornamento, ed alia 
conseguente equalizzazione adattativa, partendo dalle estremita del midambolo e 
procedendo parallelamente nei due sensi verso le estremita del burst 
preventivamente memorizzato e sincronizzato. Questo com porta u vantaggio at poier 
disporre del doppio del tempo per i calcoli. E anche possibile procedere con una 
prima passata sul semiburst di 'destra e solo dopo su quello di sinistra, usando i dati 
decisi del semiburst di destra per avere un insieme di dati decisi piu lungo su quello di 
sinistra. Questo modo di procedere migliora la precisione ma non la velocita della 
stima. 

Con gli accorgimenti posti finora in atto, e considerando che la stima congiunta 
dei parametri ottimali di equalizzazione viene estesa all'intero burst, e utile potenziare 
I'efficacia dei dati piu recenti rispetto a quelli piu lontani nel tempo, in modo da 
rendere il metodo piu reattivo rispetto alle variazioni delle caratteristiche del canale. 
Viene alio scopo introdotto un coefflciente di oblio inferiore all'unita, come d'altronde 
anche indicato nell'arte nota, che pesa opportunamente gli elementi di alcune matrici 
che memorizzano i dati storici utilizzati nella stima, prima di ogni iterazione di 
aggiornamento delle stesse. 

II grado di equalizzazione ottenuto dal metodo illustrato e generalmente 
sufficiente a consentire una corretta demodulazione nella quasi totalita delle 
situazioni trasmissive reali. Tuttavia, in presenza di interferenti non-stazionari di forte 
intensita (c / 1 <-\OdB), le decisioni errate prese durante la fase iniziale 
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deiradattamento, immediatamente dopo I'arrivo del disturbo, s&$b "Ja i 7u^ a <^ da 
superare la soglia che consente alia stima dei parametrf di equal iz^azierre ottimi di 
convergere verso valori accettabili di BER nel segnale demodulate in tempi 
ragionevolmente brevi rispetto alia iunghezza del semiburst. II forte aumento della 
probability d'errore e inevitabile a causa del ritardo intrinseco delle decisioni a breve 




termine, che rallenta I'istante di visione dell'interrerente, e queiiu ui anuazione delle 
contromisure atte a neutralizzarlo. 

Per far fronte ad interfere'nti non-stazionari di forte intensity, viene descritta una 
variante di metodo che comprende una fase di continue monitoraggio dell'energia 
ricevuta sulle antenne, e di rivelazione non coerente di eventual! picchi superiori ad 
una soglia predeterminata. Questa fase e utilmente svolta in parallelo al normale 
funzionamento. La rivelazione di un picco ad un istante n viene sfruttata per 
selezionare automaticamente il funzionamento secondo la variante. Dopodiche, per 
tutta una durata AT (tipicamente di una ventina di campioni) di un transitorio di 
convergenza delle suddette matrici dei dati storici utilizzate nella stima, vengono 
aggiornati i soli coefficienti del filtro spaziale e non i campioni stimati della risposta 
impulsiva del canale, che mantengono i valori che avevano all'istante n. Questa 
strategia d rivolta ad anticipare all'istante n la futura conoscenza sulHnterferente da 
cancellare, presente nei campioni del burst memorizzato. In tal modo llnterferente e 
"visto" in anticipo nei calcoli, e 1'algoritmo adattativo di minimizzazione dell'errore puo 
cominciare ad aggiornare i coefficienti del beamformer per cercare di eliminare il 
disturbo. Owiamente un simile anticipo non 6 pensabile per i campioni della risposta 
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impulsiva, la cui stima adattativa necessita obbligatoriamente dei dati decisi a breve 
termine. 

La stima come nella variante risulta quindi "viziata", perche per ottenere i 
vantaggi che si propone introduce una perdita di sincronismo tra i campioni dei 
segnali provenienti dalle antenne ed i relativi dati decisi. La scelta di mantenere 
cost^tlT'canipioni^U i mpu l dvanfelnra^ per lanilTcGf^ 

una misura precauzionale rivolta ad evitare propagazione di errori, altrimenti dovuti 
airaggiornamento improwiso della sola componente spaziale delle matrici strutturate 



II funzionamento originario viene automatifcamente riselezionato alio scadere 
del transitorio di durata AT . 

II vantaggio conseguito dalla variante 6 quello di rendere piu efficace la 
convergenza della stima dei parametri di equalizzazione, anche in presenza di 
improwisi e forti iriterferenti che si attivano al di fuori del midamboio, raggiungendo al 
termine di un breve transitorio una piu robusta cancellazione degli echi 



Pertanto altro oggetto d'invenzione e una variante di metodo, che attua gli 
insegnamenti appena illustrati, come anche descritto nelle relative rivendicazioni. 

E utile evidenziare che il trattamento delllnterferente come indicate nella 
variante, che necessita di avere matrici separate per la componente spaziale e 
temporale del segnale, non 6 ipotizzabile con strutture filtranti spazio-temporali come 
quelle deH'arte nota. 




di natura statistica. 




dell'interferente. 
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Volendo fare il punto sulla sequenza di fasi di equalizzazione poste in essere 



Tequalizzazione inizia in presenza della sequenza di training, una volta esaurito il 
midambolo procede utilizzando i dati stimati a basso ritardo, conviene anche 
predisporre in parallelo la possibility di attivare automaticamente il funzionamento 

fiQfinndo~" l a var i ante: Qua l ora cib dovGGGG ciwoniro; la variants ftfi fl uf j&n&- 

rapidamente i suoi effetti entro il tempo AT, alio scadere del quale conviene 
riselezionare automaticamente il funzionamento precedente, il tutto entro Tintervallo 
temporale del burst di 0,577 ms. ! 



L'equalizzatore fondato sul metodo in oggetto, e la sua variante, non 6 la 
semplice riproposizione di un dispositivo descritto in modo funzionale ricalcando i 
passi del metodo, ma bensi sfrutta le potenzialita di array sistolici realizzati mediante 
circuiti integrati digitali di tipo ASIC (Application Specific Integrated Circuit) 
appositamente configurati. 

Di conseguenza ulteriori oggetti della presente invenzione sono degli 
equalizzatori che funzionano secondo il metodo di equalizzazione oggetto 
d'invenzione, e la sua variante, e che includono degli array sistolici per la stima 
congiunta dei parametri ottimali di equalizzazione, come meglio descritto nelle 
rispettive rivendicazioni. 



Ulteriori scopi e vantaggi della presente invenzione risulteranno chiari dalla 
descrizione particolareggiata che segue di un esempio di realizzazione della stessa e 
dai disegni annessi dati a puro titolo esplicativo e non limitativo, in cui: 



dal metodo in oggetto per la demodulazione deirintero burst, risulta quanto segue: 
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in fig.1 6 indicata una schiera di antenne (array) utilizzata nella presente 
invenzione; 

nelle figg.2 e 3 sono indicati degli schemi a biocchi di altrettanti beamformer 
realizzati secondo Tarte nota; 

nelle figg.4 e 5 sono indicati dei diagram mi di radiazione sul piano azimutale, 

— ott e n^ ' ~ a PP l icand ^ 
relativamente a due differenti scenari trasmissivi in ambito radiomobile; 
in fig. 6 & indicato uno schema a biocchi di massima di un ricevitore che pone in 
particolare risalto un equalizzatore funzionante secondo il metodo di 
equalizzazione della presente invenzione; 

nelle figg.7 e 8 sono indicati due diagrammi di radiazione sui piano azimutale, 
ottenuti applicando all'array di fig. 1 un beamformer incluso nell'equalizzatore di 
fig. 6, relativamente a due scenari trasmissivi identic! a quelii che originavano i 
diagrammi di radiazione delie figg.4 e 5; 

in fig. 9 vengono mostrate due succession! temporal! di immagini affiancate, 
ordinate dall'alto verso il basso per valori crescenti dell'istante discreto n, di cui, la 
successione di sinistra indica i di diagrammi 1 di radiazione sul piano azimutale 
ottenuti applicando all'array di fig.1 il beamformer di fig. 6, e quella di destra mostra 
invece gli andamenti della risposta all'impulso all'uscita di un filtro adattato 
FADAT, incluso neirequalizzatore di fig. 6; 

nelle figg.10 e 11 vengono mostrate delle curve di BER all'uscita del ricevitore di 
fig. 6, in funzione della possibile lunghezza della sequenza di training del burst di 
trasmissione, per due diversi valori di C/l\ 
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in fig. 12 6 indicato uno schema a blocchi di massima di un ricevitore che pone in 
particolare risalto un equalizzatore funzionante secondo il metodo di 
equalizzazione della presente invenzione ed una sua variante per forti interferenti; 
in fig. 13 6 indicato un trellis relativo ad un blocco STIM incluso negli equalizzatori 
dellefigg.6 e 12;. 



in fig.t4 sono mostrati due a ndamen ti a i bbh l un g o i l bur s t , ri spe tti va m en t ^V atatatt 
alle uscite dei ricevitori delle figg.6 e 12; 




- nelle figg.15 t 16, e 17 sono indicati dei diagrammi di radiazione sul piano azM}^^^^^ 
relativi alia curva mostrata in solido in fig. 14; 

- in fig. 18 sono indicati tre andamenti di BER lungo il burst, valutati all'uscita del 
ricevitore di fig. 12, per tre diversi valori del rapporto C/l del segnale d'ingresso; 

- in fig.19 sono indicati, per confronto, due andamenti di BER lungo il burst, valutati 
aN'uscita del ricevitore di fig. 12, relativamente al metodo principale ed alia sua 
variante; 

- in fig.20 6 indicata una rappresentazione a blocchi di un blocco di stima FUNZ 
incluso negli equalizzatori delle figgi6 e 12; . 

- in fig.21 § indicata una rappresentazione a blocchi di un primo array sistolico 
incluso nel blocco FUNZ di fig.20; 

- in fig.22 6 indicata una rappresentazione a blocchi di un secondo array sistolico 
incluso nel blocco FUNZ di fig.20; 

- in fig.23 6 indicato una categoria di ulteriori blocchi inclusi nel secondo array 
sistolico di fig.22; e 



- in fig.24 sono indicati dei segnali che interessano i generici blocchi degli array 

sistolici delle figg.21 e 22. 

DESCRIZIONE DETTAGUATA 
Con riferimento alia fig.6 t notiamo lo schema a blocchi di massima di un 
ricevitore che adotta un equalizzatore funzionante secondo il metodo deirinvenzione. 
II ricevitore 6 connes su ad u n a s chiera ARY di M ant e nn a uric l orniscono d oi~ G o gnj li = 
complessivamente rappresentati mediante un vettore algebrico y(n). II vettore y(n) 
perviene rispettivamente agli ingressi di un blocco BMF ed a primi ingressi di un 
blocco FUNZ. II blocco BMF e un filtro spaziale, o beamformer, a cui perviene un 
vettore dei pesi ottimizzati w^(/z), proveniente da prime uscite del blocco FUNZ. 
Airuscita del blocco BMF sono presenti le components, reale ed immaginaria (l,Q), di 
un segnale y ouX (n) che viene filtrato da un filtro FADAT adattato alia risposta 
impulsiva del canale trasmissivo. Le componenti del segnale di banda base filtrate da 
FADAT pervengono a primi ingressi di uno stimatore STIM, alia cui uscita e presente 

una sequenza stimata {*„}. Ad un secondo ingresso del blocco FUNZ perviene una 
sequenza di training {x n } memorizzata nel ricevitore, prevent! vamente sincronizzata 
ad una copia della stessa presente nel midambolo del burst di trasmissione, 
sfruttando alio scopo metodi noti basati sul sovracampionamento del segnale ricevuto 
e la correlazione delle sequenze generate dal sovracampionamento con la sequenza 
memorizzata nel ricevitore. Da seconde uscite del blocco FUNZ fuoriesce un vettore 
di L campioni ottimizzati h opt (n) della risposta all'impulso del canale, diretti verso 
ulteriori ingressi del filtro FADAT e verso secondi ipgressi del blocco STIM. 
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II blocco FUNZ comprende a sua volta un blocco PESI che calcola un vettore di 
M pesi w(n) del' beamformer BMF, ed un blocco CAMP che calcola un vettore h(*)di 
L elementi, corrispondenti ad altrettanti campioni significativi della risposta airimpulso 
del canale trasmissivo. Comprende inoltre un sommatore ^ a due ingressi a cui 
-porx/QnnnnQ^due-segnali^numericiJrutto di opportune elaborazioni del blocchjj^ggjj^ 
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CAMP, di cui il segnale proveniente da CAMP viene preso con il segno meno. 
AH'uscita del sommatore e presente un segnale errore e n che ritorna ai suddetti 
blocchi. I primi ed i secondi ingressi del blocco FUNZ corrispondono quindi agli 
ingressi rispettivamente dei filtri PESI e CAMP, le uscite degli stessi filtri sono anche 
le prime e seconde uscite del suddetto blocco. 

La rappresentazione del blocco FUNZ evidenzia solo gli aspetti principali del 
calcolo di w^OO e h opt (n)\ una descrizione maggiormente dettagliata 

delParchitettura interna di questo blocco verra fatta illustrando la fig.20. 

Nella schematizzazione di fig. 6, per semplificare la trattazione, non sono stati 
indicati alcuni blocchi riferiti ad operazioni note compiute dal ricevitore e necessarie al 
funzionamento dell'equalizzatore, come ad esempio quelle relative alia 
sincronizzazione di portante, alia conversione/demodulazione in banda base degli M 
burst ricevuti da altrettante antenne, al campionamento ed alia conversione 
analogico/digitale degli stessi, alia memorizzazione degli M burst di simboli 
digitalizzati, al riconoscimento della sequenza di training ed alia sincronizzazione dei 
burst memorizzati, etc. Inoltre tutti gli elementi dei vettori si intendono riferiti ad 



inviluppi complessi in banda base. In particolare, il segnale y(n) proveniente dalle 

antenne e quello definito nella (21). 

Nel funzionamento, il filtro spaziale BMF e un filtro trasversale a M tappi che 
ottiene il segnale y oul (n) alia sua uscita facendo la somma pesata dei segnali y(n) 
provenienti dalle antenne e presenti nei burst memorizzati: 

II filtro adattato FADAT e lo stimatore STIM necessitano, per il loro corretto 
funzionamento, della conoscenza del vettore h„^(n), dopodiche la loro realizzazione 
e nota al tecnico del ramo. In particolare il filtro FADAT e un filtro trasversale a L tappi, 
in genere con risposta in frequenza del tipo a radice di coseno rialzato, dove L e la 
lunghezza di dispersione del canale. II blocco STIM e invece un processore che 
esegue una stima MLSE della sequenza trasmessa in base al ben noto algoritmo di 
Viterbi, applicato ad un trellis con metriche di ramo la cui espressione e quella nota 
della famiglia di ricevitori che operano con simboli pilota (sequenza di training) e 

codifica differenziale. 

Nel metodo di equalizzazione che governa il funzionamento dell'equalizzatore di 
fig.6, la stima dei vettori w^(n) e h opt (n) viene realizzata in modo congiunto, 
considerando la seguente funzione costo: 



e 2 



- ^ M L-l 
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(27) 



da cui risulta che I'errore quadratico medio e 2 , diversamente dal combinatore ottimo 
di Wiener-Hopf, risulta funzione sia dei coefficienti w i del beamformer che dei 
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campioni della risposta all'impulso h. del canale. In questo caso viene introdotta la 
dimensione temporale direttamente nel funzionale errore da minimizzare. 

Viene ora proposta una notazione matriciale, alio scopo di rendere piu compatta 
la trattazione del problema di minimo della (27). Indichiamo con Y la matrice 
corrispondente ad// campioni temporali alleM antenne: 



y T (n-i) 




y T (n-N + l) 

! 

ottenuta dal vettore y(*)- [ yi (n\ y 2 (n),....,y M (n)J i cui M element! sono i campioni 



dei segnali ricevuti sulle antenne; per cui: 



y„.N+i y n -N+i 



y n+ M-2 



(29) 



Per quanto riguarda il numero N, si vedra quando verranno illustrate le figg.10 e 
1 1 che un valore compreso tra 20 e 30 (N = 26 per il midambolo GSM) 6 sufficiente a 
fornire una buona base statistica per la stima dei parametri di equalizzazione ottimi. 

Indichiamo con X la matrice dei dati trasmessi (per ora limitati alia sola 
sequenza di training), ciascunp con duratai: 



x r (") 
x T («-l) 

x T (n-N + l) 



(30) 



ottenuta dal vettore \(n)^[x n (n), x n+l (n\....,x^ L _ l (n)J i cui L elementi sono dei 
dati trasmessi in corrispondenza della sequenza di training {*„}, e opportunamente 



m 



emorizzati, assumendo una lunghezza di dispersione del canale pari aL; per cui: 



-X— 



n+l 



e con: 



X n~N+l X n-N+2 '" X n-N+L 



(32) 



^[^,....,^7 ( 33 ) 
sono rispettivamente indicati i vettori dei coefficienti del filtro spaziale BMF e della 
risposta all'impulso del canale. 

II problema proposto neil'equazione (27) puo essere quindi riscritto in forma 



matriciale, come segue: 

(w^,h op ,) - argmine 2 (w,h) = arg min||Yw - Xh|| 



imp 



(34) 



(35) 



dove nella (35) viene introdotto il vincolo sull'energia della risposta all'impulso per 
evitare la soluzione banale (w opt - 0, h opt = o) . II vincolo (35) e una conseguenza del 
fatto che i bit della sequenza trasmessa sono assunti come temporalmente incorrelati, 
valendo altre ipotesi la (35) dovrebbe essere modificata in modo congruente. 
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II problema del minimo vincolato pu6 essere affrontato applicando il metodo di 
Lagrange. Alio scopo, detto y un parametro reale (moltiplicatore di Lagrange), 
consideriamo il seguente funzionale: 

" j(w, h, X') - || Yw - Xh|| 2 + y (IMI 2 - l) (36) 
La ricerca del minimo vincolato nella (34) si riduce quindi alia ricerca di una 



"estremante" (w, h) del problema di minimo libero del funzionale j(w,h,y). La 
suddetta estremante pu6 essere ricercata considerando il sistema: 

^ - / w (w,h, Y ) - 2Y'(Yw-Xh) - o"' (37) 

|£ = / h (w,h, Y ) - -2X'(,Yw-Xh)+2 Y h ~ 0 (38) 
on 

|£ = J x (w,h, Y ) - h'h-1 = 0 (39) 

L'equazione (37) pud essere interpretata come I'introduzione del principio di 
ortogonalite dell'errore, 6 infatti possibile esprimere l'equazione (37) come: 

Riprendendo, per comodit&, alcune espressioni in parte gia introdotte, 
indichiamo con la matrice di covarianza spaziale espressa da: 



R„ = £[y(0y*(r)] (41) 
ed analogamente con e R , rispettivamente, le matrici di covarianza sui dati e 



di cross-covarianza tra dati ed ingressi alle antenne espresse da: 



R„ - E[x(t)x\t)] (42) 
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- E[x(t)y'(t)] (43) 
Nel seguito useremo le matrici di covarianza campionaria definite come segue: 

1 " 




) 



iV n-l 



(17) 



Possiamo pertanto definire una matrice R strutturata come segue: 



(45) 



R 



R. 



yy 



R xy 

R„ 



(46) 



che ci consentira una piu semplice rappresentazione del problema. Le equazioni (37) 
e (38) possono essere rispettivamente riscritte come: 

Rw = R\h (47) 



R v w-R„h = yh 



(48) 



Dall'equazione (47) e possibile ricavare w in funzione di h e sost'rtuirlo nella 
(27) per ottenere I'errore quadratico medio. Sostituendo invece il vaiore di w nella 
(48) e possibile dimostrare la seguente uguaglianza: 

e 2 (w,h)=h'R, i h = -y (49) 

dove: 



(50) 



6 detto complemento di Schur di in R. La stima del vettore h(#i) che minim izza 



I'errore 



quadratico medio, corrisponde quindi alia stimai dell'autovettore h opi (») della 



matrice associato all'autovalore y a modulo minimo. 

E possibile quindi stimare anche i coefficienti w opr («) del beamformer 




Riassumendo, la scelta dei pesi ottimi che minim izzano la 
essere ricondotta a questi punti fondamentali: 

1. calcolo delle matrici R w , e del complemento di Schur R^ ; 

2. stima dell'autovettore della matrice R^ associato airautovalore a modulo minimo 

3. stima dei coefficienti w opr (/i) del beamformer. 

Assumendo di realizzare il metodo di equalizzazione riassunto ai punti 1., 2., 3., 
mediante un programma che controlla un elaboratore elettronico a precisione finita, il 
calcolo del complemento di Schur pud risultare un'operazione assai critica, 
soprattutto se la lunghezza di parola dell'elaboratore e molto limitata. 

La suddetta criticita viene notevolmente ridotta grazie alPutilizzo dei metodi 
sviluppati in una teoria sulla decomposizione parziale in forma triangolare delle 
matrici, nota in letteratura con la denominazione : di "Displacement Theory", o anche 
•Displacement Structure", riportata nella bibliografia in seguito indicata, a cui si 



rimanda. Elemento chiave di questa teoria e un'entita denominata "fattore di 
Cholesky" di una matrice, nel caso in oggetto la matrice R, il cui calcolo e svolto 
mediante il semplice aggiornamento di un valore precederrte, intervenendo con 
opportune trasformazioni (rotazioni di Givens) direttamente sui nuovi dati. Come si 
vedra. raggiornamento del "fattore di Cholesky" fornisce una matrice, strutturalmente 



diversa dalla R ma con il medesimo contenuto informativo, i cui elementi facilitano 
grandemente il calcolo di w opt e h opt . Viene in tal modo escluso il calcolo esplicito 
9^ delle matrici di covarianza campionarie ed il ricorsp ad operazioni critiche, quali quella 

dell'inversione matriciale diretta, poiche le suddette rotazioni non comportano alcun 
calcolo del determinants Questo modo di agire comporta un migliore 
condizionamento del problems di calcolo, ed una conseguente riduzione della 
lunghezza minima di parola deil'elaboratore. 

Un tale approccio consente I'ulteriore vantaggio di poter estendere il calcolo 
all'intero burst; la soluzione descritta e infatti intrinsecamente ricorsiva (iterativa) . Per 
contro, un metodo non ricorsivo, che comporti cioe il calcolo esplicito delle matrici di 
covarianza campionarie, risulta caratterizzato da un'elevata complessita nel calcolo 
del complemento di Schur R, 1 mediante la (50), che rende difficoltosa I'estensione di 

un tale metodo all'intero burst. 

Occorre sottolineare che gli esempi applicativi della "Displacement Theory" 
finora noti, riguardano esclusivamente I'ottimizzazione delle prestazioni di 
equalizzatori solo temporali. Elemento di novita e stato quello di ricondurre la 
formulazione matematica della soluzione del problema di minimo del funzionale errore 
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utilizzato nella presente invenzione, alia formulazione fornita da tale teoria per la 
soluzione del problema di minimo in equalizzatori MMSE-DFE solo temporali. In virtu 
dell'analogia evidenziata tra le due formulazioni e stato quindi possibile sfruttare tutta 
la conoscenza sviluppata per gli equalizzatori temporali adattandola aU'equalizzatore 
spaziale e temporale dell 'invenzione. 



Cio posto, consideriamo la matrice R descritta nella (46) la quale, essendo 
definita positiva, garantisce I'esistenza di un unico fattore di Cholesky L tale che: 

R - LL* (52) 

Uno dei metodi fondamentali per la decomposizione in fattori triangolari, e il 
processo di riduzione di Schur descritto: 

• nell'articolo di J. Chun e T. Kailath, intitolato "Generalized Displacement Structure 
for Block-Toeplitz, Toeplitz-Block, and Toeplitz-deived Matrices", pubblicato in 
NATO ASI Series, vol. F 70; 

• nel volume intitolato "Numerical Linear Algebra Digital Signal Processing and 
Parallel Algorithms", autori G. H. Golub e P Van Dooren, edito da Springer-Verlag, 
1991. 

II processo consente, dopo M passi di decomposizione, la riduzione della 
matrice R nella forma: 



R 



[V If]- 



0 0 

0 R, x 



(53) 



dove R^ pud essere quindi considerato come il residuo della decomposizione 
parziale in forma triangolare. Questo approccio ci consente di calcolare la matrice 



V:: 
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senza dover invertire la matrice di covarianza spaziale R^. Partendo dal 
complemento di Schur, derivato dal processo di decomposizione, 6 possibile stimare 
h con metodi classici di stima deirautovettore associato all'autovalore a modulo 



*opt 



minimo, ad esempio, il metodo detto delle potenze inverse. 



che consentono la riduzione dell'equazione (51) nella forma: * 

L'w^ = U'h^ (56) 
La forma triangolare inferiore della matrice L, come vedremo in seguito, 
consente la soluzione dell'equazione (56) con un semplice processo detto di 
"back-substitution". 

Analizziamo ora i diversi passi sequenziali del metodo di equalizzazione che 
utilizza nei calcoli il fattore di Cholesky, riportando per ognuno la relativa complessita. 

• passo 1: applicazione di un algoritmo noto di riduzione veloce di Schur della 
matrice R per la decomposizione in fattori di Cholesky; 

♦ passo 2: applicazione di un algoritmo noto ! in letteratura come "Metodo delle 
potenze inverse" per la stima deirautovettore (h op( ) di R, 1 associato all'autovalore 
a modulo minimo; 



Italtel 




decomposte le matrici R w e R v nella loro forma triangolare: 

: R w = LL* (54) 
R = UL' (55) 



; 3 



@ 
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passo 3: risoluzione. dell'equazione (56) per la stima dei coefficienti ottimi del 
beam-former (w opr ), mediante un metodo noto in letteratura come 
"back-substitution" 



Passo 1 

n^iaJeoria-deJle,ro.atricLdj Toeplitz e noto che una generica^^nce 




una strutt'ura spostamento tipo-Toeplitz se la matrice differenza VR(«) definita come: 

VR(n) = R(n) - ¥(n)K(n - l)F' («) (57) 
ha basso rango (rank(VR) « rank(R)). La matrice R espressa nell'equazione (46) 
rientra in questa categoria e, in accordo alia suddetta teoria, pu6 essere calcolata in 
mode iterativo, ponendo arbitrariamente R(o)-0, ed estendendo a tutto il 
midambolo I'aggiornamento definito come segue: 

R(«)- R(^-l)^g n g; (58) 
dove il vettore aggiornamento g n , definito come: 



Sn 



(59) 



6 detto generatore della matrice differenza, in qupnto imponendo F(n) = I possiamo 

scrivere VR(/i) = g n g* n . 

Alia stessa stregua § possibile aggiornare il fattore di Cholesky della matrice 
R(/x), grazie alia speciale natura tempo-variante della stessa. Detto quindi L(/i) il 
fattore di Cholesky della matrice R(n), tale che R(n)-L(/i)X7(/i), e possibile 
verificare che: 
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L(//)l/(») = [L(»-0 g„] 



r(/;-i) 



(60) 



Questo suggerisce un metodo efficiente per raggiornamento del fattore di 
Cholesky attraverso I'uso di trasformazioni unitarie, di modo che, partendo dal fattore 
calcolato all'istante n-1 e dal generatore all'istante n, sia possibile determinare il 



.fattore all'istante n senza dover caicoiare la maince R{n) . 



nuovo 



E infatti possibile, posto arbitrariamente L(0) = 0 , esprimere ii fattore 
aggiornato come: 



[E(«) o] = [l(»-i) g„]Q(») 



(61) 



dove Q(«) e una rotazione unitaria, in genere attuata come una sequenza di 
trasformazioni elementari, Q 0 ,Q,,-", che, applicata a [l(«-1) g„] produce 
I'annullamento della colonna del generatore g„ introducendp uno zero per ogni 
trasformazione elementare. 

Come esempio di aggiornamento ricorsivo, siano / 0 («) e/ 0 (»-l) i primi 

elementi in posizione (o,o) rispettivamente delle matrici L(«) e L(« - 1) , e g„ 0 il 
primo elemento del vettore del generatore. Ne deriva che la prima riga della matrice 



L(w) puo essere scritta come: 

[/ 0 W 0 0] = [/„(«- 1) 0 g„; 0 ]Q 0 
Se applichiamo la stessa trasformazione all'intera matrice ne deriva: 



(62) 



[L(«-l) g„]Q 0 = 



hi") 



L,(//-l) g„., 



(63) 




Itatte! 



dove L,(/7-l) e g nl sono le matrici ottenute dopo aver eliminato il primo elemento 



della colonna del generatore. 



Se procediamo ricorsivamente neiraggiornamento deU'intero fattore L(/7) , lo 



possiamo vedere decomposto nella seguente forma triangolare: 



'Un) 0 



■am 



■U(/i) L» 



dove L e U sono le matrici definite nella (53), e L, e il fattore di Cholesky di R, 1 , 
tale per cui: L,L* = R s x Come si pud notare dalla (64), le matrici L e L, sono 
triangolari mentre la U e necessariamente rettangolare, o eventualmente quadrata. 

La fattorizzazione della (64) risulta fondamentale in quanto consente il calcolo 
del vettore dei pesi w 0 ^(/?) senza dover effettuare flnversione della matrice 
(vedi equazione 50) e, come vedremo tra breve, senza dover invertire la matrice , 
operazione questa richiesta dall'algoritmo delle potenze inverse per il calcolo 
deirautovettore ( h 0pt ) associato airautovalore minimo. 

La decomposizione ricorsiva per I'aggiornamento deirintero fattore di Cholesky 
richiede un numero di moltiplicazioni complesse pari a 2(M + L)(M + L - 1) ed un 
numero di somme pari a (M + L)(M + L - 1). 



II metodo delle potenze inverse per il calcolo dell'autovettore associato 
airautovalore minimo richiede la soluzione in modo iterativo dell'equazione: 



Passo 2 



(65) 



• 



♦ 
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KM = ^ 

dove la normalizzazione della (66) viene introdotta come risultato del vincolo 
suirenergia della risposta all'impulso del canale stimato: |h| 2 = 1. A causa del 
suddetto vincolo conviene scegliere un vettore iniziale h 0/ ,,(0)con almeno un 
elemento non nullo. 

La riduzione in fattori di Cholesky consente il calcolo del suddetto autovettore 
senza dover ricorrere all'inversione della matrice, in quanto e possibile riscrivere la 
(65) nella forma: 

L,(/7)L>)h n = h 0/ >-l) (67) 
L'equazione (67) pud essere scomposta in due fasi elementari, introducendo un 
vettore intermedio u n definito come: 

L»h„ = Un (68) 
in modo da poter esprimere u n come soluzione di: 

L,(»K = h opt (ti-\) (69) 
Possiamo risolvere I'equazione (69) rispetto a u n e quindi, avendo trovato u rt , 
siamo in condizioni di risolvere I'equazione (68) rispetto a h„. Entrambe le equazioni 
sono semplici da risolvere in quanto le matrici interessate sono matrici triangolari. Se 
dovessimo risolvere equazioni con generiche matrici quadrate LxL, ci dovremmo 
aspettare un numero di operazioni aritmetiche deirordine di L 3 . Essendo invece le 
matrici triangolari e possibile raggiungere la soluzione in sole L 2 operazioni. Per 
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illustrare tale approccio, noto in letteratura e detto di back-substitution, sviluppiamo le 



# 



equazioni (69) e (68): 




l U h r,.L = U L 
1 Il-Kl + l l-\.L-A.l-\ = 



Partendo dal primo termine della (70) possiamo calcolare «, : 

e sostituirlo nelle altre L - 1 equazioni. Questo mostra che possiamo ridurre il 
problema ad un sistema triangolare di dimensioni (L - 1) x (L - 1) ad un costo di L - 1 
operazioni di moltiplica-e-somma ed una divisione. Con questo metodo possiamo 
risolvere il sistema al costo di L divisioni e (L - 1) + (L - 2) + ... +1 operazioni di 
moltiplica-e-somma, per un totale quindi di VzL{L - 1) operazioni di moltiplica-e-somma 
e L divisioni. Dopo aver risolto I'equazione (69) rispetto a u n possiamo risolvere 
Tequazione (68). II numero di operazioni richiesto e del tutto identico a quello 
necessario per la soluzione del primo insieme triangolare di equazioni. 



Passo 3 
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Anche i! calcolo dei coefficienti ottimi del beamformer puo essere effettuato per 
back-substitution nella (56), direttamente nello spazio dei segnali (senza calcolo 
esplicito della matrice di covarianza), calcolando: 



V(n)w opt (n) = V\n)h opt (n) 



(56) 



In questo caso il precedente passo 2 deve essere esteso introducendo anche il 
calcolo del secondo termine della (56), per un numero totale di ML + !4M2 operazioni 

di moltiplica-e-somma e M divisioni. 

La seguente tabella riassume i dati relativi alia complessita di calcolo nei tre 
passi illustrati, riportando anche la potenza di calcolo di un ipotetico elaboratore 
(espressa in MIPS) richiesta per I'esecuzione in tempo reale. Viene considerato un 
aggiornamento dei coefficienti, esteso al solo midambolo, per una schiera di antenne 
di M = 8 ricevitori e per un filtro adattato con lunghezza L = 5. 



Passo n° 


N° moltiplicazioni complesse 


N° somme complesse 


MIPS 


1 


2(M + L){M + L- 1) 


(M + L)(M + L- 1) 


8,97 


2 


L(L-1) + 2L 


L(L-1) 


0,86 


3 


ML + 'AM 2 


ML + ViM 2 


0,098 



Le potenzialita di una stima congiunta possono essere dimostrate considerando 
lo scenario precedentemente proposto in fig. 4, dove la dispersione temporale 
introdotta dal canale e rappresentata da un eco del segnale ritardato di 5 jus. 

La fig.7 mostra il diagramma di radiazione sul piano azimutale determinato dal 
beamformer BMF di fig.6, dove risulta evidente il diverso comportamento rispetto alia 
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situazione di fig. 4 nei confronti dell'eco ritardato. II diagramma di radiazione di fig. 7 
presenta infatti un livedo di energia nella direzione dell'eco, della stessa entita del 
livello nella direzione dell'utile. Cio consente, a prezzo di un'interferenza 
intersimbolica controllata, di mantenere alta I'energia del segnale utile presente nel 

processo di stima della sequenza trasmessa [x n ] . 



La robustezza del nuovo funzionale proposto pud essere dimostrata in fig. 8, 
proprio considerando il comportamento del beamformer BMF a fronte dello stesso 
scenario critico di fig. 5, in cui era presente una replica deirinterferente cocanale nella 
stessa direzione del segnale utile. 

Come si puo notare dal diagramma di radiazione di fig. 8, I'energia e ora 
trasferita prevalentemente nella direzione dell'eco dell'utile, cosicche il diagramma 
presenta uno zero anche nella direzione dell'utile, oltre che dell'interferente. Questo 
comportamento e spiegabile in quanto I'introduzione della dimensione temporale nel 
funzionale consente la combinazione degli echi, anche qualora questi dovessero 
essere temporalmente ritardati. Possiamo notare la capacita del funzionale misto di 
selezionare nello spazio, tra i diversi echi ritardati che sono coerenti con la sequenza 
di training, quello caratterizzato dal maggiore rapporto tra segnale ed interferente. 

Questo comportamento consente di gestire anche utenti mobili con direzioni 
nominali sovrapposte, qualora abbiano una opportuna dispersione spaziale degli echi. 
Si puo dimostrare sperimentalmente come, in ambito RCS, cio consenta una notevole 
riduzione della probability di handover intracella. 

La stima congiunta dei parametri di equalizzazione ottimi, necessita tuttavia di 
un maggior numero di parametri rispetto ad una stima solo temporale. In particolare, 
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nel caso dei sistemi FDMA-TDMA con sequenze di training relativamente corte (come 
appunto nel GSM), si presentano due potenziali problemi: 

• una crescita della complessita di calcolo nel metodo di equalizzazione congiunta in 
oggetto, rispetto ad un'equalizzazione disgiunta dell'arte nota; 

• una maggiore varianza nella stima dei parametri di equalizzazione ottimi, e quindi 
una piu lenta velocita di convergenza della stessa. 

La lunghezza della sequenza di training rappresenta un parametro critico negli 
equalizzatori utilizzati, o utilizzabili, nei sistemi di cui sopra, in quanto gli algoritmi per 
la stima congiunta dei parametri ottimali di equalizzazione devono identificare la 
risposta all'impulso.del canale, filtrata dal beamformer BMF, disponendo di un numero 
di campioni limitato, ovvero i 26 del midambolo del burst GSM. 

La fig. 9 mostra, partendo dall'alto verso il basso, una sequenza di diagrammi di 
radiazione sul piano azimutale e di rappresentazioni della risposta all'impulso del 
canale, "fotografate" al compimento di n iterazioni della loro stima congiunta, compiuta 
mediante gli algoritmi adottati ai passi 1, 2 e 3. Lo scenario raffigurato e relativo ad un 
segnale utile e ad un interferente cocanale (indicati rispettivamente con "o" e V\ 
sull'arco piu esterno), entrambi accompagnati da un proprio eco ritardato di 10 
(indicati con gli stessi simboli posti appena al di sotto deirarco piu esterno). 

Con riferimento alia fig. 9 e possibile notare come la convergenza complessiva 
sia dominata dalla velocita di convergenza della risposta all'impulso h opt (n)- 
ANIstante n = 6, quando la base statistica non e sufficientemente estesa per 
effettuare una corretta stima della risposta all'impulso, il beamformer BMF introduce 



quando cioe la stima di h opl (n) porta a riconoscere il solo 



degli zeri anche in com'spondenza degli echi del segnale 
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beamformer BMF non si comporta ancora in modo ottimo, in quanto introduce uno 
zero in posizione corrispondente all'eco ritardato del segnale. Possiamo pero notare 
che gia per n = 26, cosi come per n = 30, I'identificazione del canale consente 
I'individuazione dell'eco ritardato come segnale utile, e quindi il diagramma di 
radiazione presenta degli zeri solo in com'spondenza dell'interferente e del relativo 



Una piu dettagliata analisi della velocita di convergenza pud essere svolta 
valutando il BER all'uscita del ricevitore di fig.6 al variare della lunghezza del 
midambolo, quando esso opera in condizioni critiche, ad esempio con C/l = - 10 dB. 

Le figg.10 e 11 mostrano appunto i risultati di questa analisi riferiti ad uno 
scenario che comprende un interferente cocanale, oltre che all'utile, per tre diversi 
valori del rapporto segnale-rumore SNR all'ingresso del ricevitore: 10, 20, e 30 dB. 
Piu precisamente, le curve hanno in ascissa la lunghezza della sequenza di training, 
espressa come numero di campioni, ed in ordinata il BER. La fig. 10 e riferita ad una 
situazione cbn C/l = - 10 dB, mehtre la fig.1 1 ad un C/l = 0 dB. 

Con riferimento alle figg.10 e 11, si pud notare che in ambedue i casi le 
prestazioni del ricevitore si possono considerare stabilizzate per valori di lunghezza 
del midambolo compresi tra 20 e 26 campioni. In questa zona e infatti localizzata la 
quasi totale variazione di pendenza delle curve, il che indica la raggiunta convergenza 
della stima dei parametri di equalizzazione ottimi all'interno della sequenza di training. 



eco. 
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Con riferimento alia fig. 12, notiamo lo schema a blocchi di massima di un 
ricevitore il cui equalizzatore e quello di fig. 6 con I'aggiunta di una parte che tiene 
conto del funzionamento adattativo fuori midambolo. 

Come si pud notare confrontando le figg.6 e 12, le due rappresentazioni 
deU'equalizzatore differiscono per il fatto che in fig. 12: 



il: 



- e evidenziato un blocco aggiuntivo RIT che introduce un ritardo D sul percorso del 
bus del vettore d'ingresso y(//) verso il blocco FUNZ; 



- al 



blocco FUNZ perviene un vettore x'(/?) invece della sequenza di training {*„}; 



- c'e un nuovo blocco DECTEMP che memorizza il vettore x'(n), uscente dal blocco 
STIM. 

Come si vedra illustrandone il funzionamento, la schemartizzazione di fig. 12 e 
idonea al funzionamento in presenza di canali con caratteristiche tempo-varianti e/o in 
presenza di interferenti non-stazionari, che si attivano cioe al di fuori della sequenza 
di training, e che quindi non possono essere individuati e cancellati considerando il 
solo midambolo nella stima dei parametri ottimi deirequalizzatore. Proprio per questo 
le due schematizzazioni delle figure 6 e 12 sono tra loro complementary nel senso 
che durante il processamento del midambolo e attiva la prima, ed al termine, e per 
tutta la durata rimanente del burst e attiva la seconda. In pratica ambedue le 
schematizzazioni sono riferibili a due diverse modalita operative di un'unica struttura 
di equalizzatore. 

Come e noto, I'effetto Doppler e una delle cause che in un contesto radiomobile 
contribuisce maggiormente a rendere tempo-varianti le caratteristiche dei canali. Esso 
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consiste in una piccola variazione di frequenza del segnale conseguente al 
movimento del mobile, che comporta una lenta variazione temporale dell'ampiezza 
degli echi ricevuti. II canale risulta cosi tempo-variante. e anche se le variazioni 
avvengono molto lentamente, puo essere utile cercare di inseguirle aggiornando man 
mano i valori ottimali di parametri h opt e w opl stimati durante il midambolo. In tal 




modo..si riesce a migliorare la qualita di demodulazione, poiche lo stimatore STIM e 
reso adattativo, in quanto ricalcola ad ogni passo le metriche di ramo in base alia 
nuova stima di h opl , che a sua vola risente di quella di w 

opt ' 

Operativamente. puo risultare conveniente aggiornare (nel modo che si vedra) 
la coppia ottimale di parametri h opl e stimati durante il midambolo, ed in 
concomitanza equalizzare il segnale ricevuto. partendo dalle due estremita del 
midambolo e procedendo parallelamente nei due sensi verso le estremita del burst, 
preventivamente memorizzato e sincronizzato. II maggior parallelismo circuitale che 
ne deriva consente di avere a disposizione un tempo doppio per I'esecuzione dei 
calcoli. Questo e un vantaggio non da poco. visto I'alto numero di singoli element! di 
calcolo da eseguire in un tempo di simbolo di soli di 3,69 us, che verrebbe 
raddoppiato utilizzando due semiburst in parallelo. 

La presenza di interferenti non-stazionari all'interno del bus di trasmissione e 
invece spiegabile con il fatto che le BS relative a celle diverse non sono 
necessanamente sincronizzate fra loro. Nel caso, ad esempio, di un approccio RCS 
al problema dell'aumento dell'efficienza spettrale. ci si puo quindi aspettare di avere 
interferenti cocanale in celle molto vicine a quella in cui si trova il mobile che trasmette 




il segnale utile. Poiche non e possibiie sincronizzare tutte le BS, e dunque probabile 
che un interferente disturbi la trasmissione solo durante una parte del burst e che 
quindi non sia "identificabile" e cancellable durante la trasmissione della sequenza di 
training. 

Lequalizzatore di fig.12, grazie alia possibility di aggiornare durante I'intero 
burst fa coppia di parametri ottimali h opt e w opt stimati durante il midambolo, e posto 
cosi in grado di affrontare anche il caso in cui uno o piu interferenti disturbino il 
segnale utile a partire da un istante esterno alia finestra temporale in cui viene 
trasmessa la sequenza di training. 

Nel funzionamento, si pud gia fin da ora anticipare che, con le precisazioni che 
seguiranno, I'algoritmo svolto dal blocco FUNZ rimane quasi del tutto invariato. Cio e 
anche intuibile dalla grande similitudine esistente tra le raffigurazioni delle figg.6 e 12, 
resa ancora piu evidente dal fatto che il blocco RIT che ritarda di un tempo discreto 
D il vettore y(n) pud essere incluso nel blocco FUNZ. 

II funzionamento adattativo richiede la conoscenza della sequenza dei dati 
trasmessi {x„}. Quando questi non sono noti, ovvero quando non e piu disponibile la 

sequenza di training, occorre sostituirii con i dati {x n } decisi nel processo di 
demodulazione. Viene in tal modo introdotto un inevitabile ritardo tra i campioni dei 
segnali alle antenne ed i corrispondenti dati della sequenza stimata MLSE, necessari 
entrambi per determinare Terrore. Infatti, come si pud notare dalla fig.12, prima che 
sia disponibile il dato x n relative all'istante n , bisogna attendere sia il ritardo L 
introdotto dal filtro adattato FADAT, sia il ritardo dovuto al tempo di latenza 
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dell'algoritmo di Viterbi necessario aU'emissione di dedsionj^f^T^^^^ :dabi!i. 

I'll 

Piu precisamente, se D e il ritardo complessivo (L + latei\za di ^terbij^i^tl decisi 
della sequenza {x„} sono disponibili solo all'istante w + Z), e quindi la stima dei 
parametri di equalizzazione ottimi h opt (n) e w opt (n) all'istante ti pud sfruttare la 
- con nr.nfinzaHjeqli-in^ 



ovvero i vettori y(n - D) e x(n-D). II blocco RIT compensa il ritardo gia esistente 

sull'altra via. - 

E utile ribadire che il vettore x(w) contiene ora dei dati decisi x n , mentre durante 
la stima congiunta dei parametri ottimali di equalizzazione, fatta durante il midambolo, 
lo stesso vettore conteneva dati gia noti x n e quindi privi di ritardo. L'affermazione 
fatta sul vettore y(n- D) significa che, per quanto concerne il calcolo di h opt (n) e 
w in) la media all'istante n -esimo calcolata sugli N elementi delle colonne delle 

opt v / » 

matrici R^, R^, e R^ riguarda elementi compresi tra I'istante n-D e I'istante 
n-(D + N). 

Cio che in pratica si verifica a causa del ritardo D, e che la stima congiunta dei 
parametri ottimali di equalizzazione all'istante n risulta in ritardo rispetto alle variazioni 
del canale, che si ripercuotono senza ritardo alcuno sui simboli y(n) ricevuti. Si 
scontrano quindi due necessita. La prima e quella di mantenere alto il tempo di 
latenza dell'algoritmo di Viterbi, per garantire decisioni corrette che evitino 
propagazione di errori nell'anello di stima. La seconda e quella di fornire i dati decisi 
nel piu breve tempo possibile al blocco FUNZ, per non introduce un ritardo eccessivo 




neiraggiornamento dei parametri ottimali di equalizzazione, tale che pregiudichi la 
possibility d'inseguire variazioni veloci nelle caratteristiche del canale. 

L'approccio seguito nel metodo in oggetto e quello di utilizzare nella stima dei 
parametri ottimali di equalizzazione, una sequenza di dati stimati a basso ritardo, 
rappresentata dagli L elementi del vettore x'(//). Questi e ottenuto forzando lo 
stimatbre STIM a decidere gia dopo pochi passi (circa 4 o 5), qual e la sequenza piu 
probabile, ottenuta risalendo nel trellis il percorso a massima verosimiglianza. I dati 
definitivamente decisi della sequenza {x(n)} si rendono invece disponibili alia fine del 
burst, in modo che siano il piu possibile affidabili. 

La tecnica di cui sopra, che risolve il problema di poter disporre 
preventivamente dei dati di una sequenza stimata quasi-correttamente, e detta di 
"Per-Survivor Processing", o PSP, e viene descritta nell'articolo intitolato "Per-survivor 
processing: a general approach to MLSE in uncertain environments", autori: R. 
Raheli, A. Polydoros e C. K. Tzou, pubblicato in IEEE Trans. Commun., vol. 43, 
pp.354-364, February/March/April 1995. Mediante questa tecnica si pud utilizzare 
I'algoritmo di Viterbi in presenza di uno o piu, parametri ignoti, o non ben determinati, 
in questo caso i campioni della risposta impulsiva ottimizzata, che se fossero invece 
noti potrebbero essere vantaggiosamente utilizzati neiraggiornamento delle metriche 
di ramo. 

Poiche le decisioni a basso ritardo possono risultare errate, si puo generare un 
fenomeno di propagazione degli errori dovuto ad un circolo vizioso che si instaura 
nell'anello formato dalla retroazione di x'O?) verso il blocco FUNZ. Considerando 
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infatti gli elementi del vettore x*(/?) come il frutto di singole decisioni a basso ritardo a 
cadenza di simbolo, come appartenenti ad una sequenza che scorre in un registro a 
scorrimento lungo L (questo e quello che effettivamente avviene nella 
rappresentazione degli stati secondo Viterbi), una decisione errata finirebbe per 
produrre i suoi effetti negativi per un numero di tempi di simbolo uguale a quello 



stabilito per le decisioni a basso ritardo. 

II metodo di equalizzazione in oggetto elimina anche questo secondo 
inconveniente sostituendo in blocco, ad ogni tempo di simbolo /?, tutti gli L elementi 
del vettore x'(/?) con degli altrettanti dati temporanei ricopiati direttamente dal 
percorso a massima verosimiglianza nel trellis valido al tempo di simbolo /?-esimo. 

La fig. 13 mostra appunto un esempio in cui all'istante N viene presa nel trellis 
una decisione temporanea sbagliata, che pero non si propaga perche all'istante 
successivo JV+ 1 la sequenza decisa che viene ricopiata e quella corretta. Non 
sempre in realta le cose possono andare in tal senso, ma con buona probability cio e 
verificato in presenza di interferenti non particolarmente intensi. 

La modalita operativa appena illustrata rende quindi piu accurato 
Taggiornamento dei parametri ottimi al di fuori del midambolo, permanendo 
I'approssimazione introdotta dalle decisioni a basso ritardo. Inoltre, per migliorare il 
comportamento adattativo dell'equalizzatore di fig. 12, occorre introdurre un 
coefficiente A cosiddetto di "oblio", in modo da pesare meno i dati meno recenti. Con 
cio, ricordando la (58), la matrice R all'istante // pud essere scritta come: 

R(//) = AR(//-l) + g„ g ; (72) 
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y(7 ° (73) 



x'(n) 

dove A e una costante reale arbitraria tale per cui: 0 < X < 1 , e g„ e un vettore di 
aggiornamento definito come nella (59), con la differenza che i suoi elementi sono qui 
le decisioni a basso ritardo. 

Alio stesso modo, ricordando la matrice differenza (57) e la (58), si pud scrivere: 

R(„) = F(w)R(«-l)F-(/j) + g„g; (74) 

essendo: 

F(«) = VTl (75) 
Con questa posizione I'aggiornamento del fattore di Cholesky pud essere 
derivato direttamente dall'equazione (74). Detto infatti L il fattore di Cholesky di R , 
possiamo aggiomarlo, analogamente a quanto fatto nella (61), seguendo la relazione: 

0] = yiL(n - 1) g.]Q(») (76) 

Con I'insieme delle precisazioni di cui sopra, la stima dei parametri h gpl (n) e 

v/ opt (n) , resa adattativa, e ottenuta applicando I'algoritmo ricorsivo gia illustrato per la 

stima non adattativa compiuta dal blocco FUNZ di fig. 6. 

Facendo riferimento alle figure 14, 15, 16, e 17, vengono dimostrate le effettive 
capacita di equalizzazione adattativa del metodo in oggetto. 

La fig. 14 mostra due andamenti di BER lungo il burst di altrettanti segnali 
demodulati, per un interferente che disturba fino al 50° simbolo. Un primo andamento 
(a tratteggio) e relativo al metodo di equalizzazione in oggetto, supponendo di 
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congelare la stima al di fuori del midambolo, mentre il secondo andamento (in solido) 
e relativo alio stesso metodo consentendone I'esplicazione del funzionamento 
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adattativo. II segnale d'ingresso presenta —=\0dB e —-OdB\ gli andamenti 



sono mediati su 500 realizzazioni. II BER indica la probability che il generico 




notare dalla prima curva, se si demodula il burst utilizzando i parar 
equalizzazione stimati con la' sequenza di training, le prestazioni degradano 
notevolmente per tutto il tempo in cui e presente I'interferente. La seconda curva 
mostra che il funzionamento adattativo, invece, riesce a limitare il degrado delle 
prestazioni. Infatti, procedendo dal centro verso le estremita, si puo vedere come le 
prestazioni degradirio in corrispondenza dell'arrivo (nel calcolo) dell'interferente, per 
poi migliorare in seguito all'equalizzazione. Cio significa che I'equalizzatore e riuscito 
ad adattarsi alia nuova situazione trasmissiva, ed a variare il diagramma di radiazione 
della schiera di antenne in modo da cancellare I'interferente. Le successive figg.15, 
16, e 17 mostrano dei diagrammi di radiazione sul piano azimutale relativi alia curva in 
solido di fig. 14, 

La fig. 15 mostra il diagramma di radiazione stimato alia fine del midambolo, 
quando Tinterferente non e ancora presente. Le direzioni di arrivo degli echi del 
segnale utile sono. rappresentate in solido con tratti di linea di lunghezza 
proportionate alia potenza del rispettivo eco, espressa in dB. 
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La fig. 16 mostra I'evoluzione che subisce il precedence diagramma di radiazione 
subito dopo Tindividuazione dell'interferente, i cui echi sono rappresentati da linee 
tratteggiate di lunghezza proporzionale alia potenza del rispettivo eco, espressa in dB. 

La fig. 17 mostra il diagramma di radiazione quando la stima congiunta dei 
parametri di equalizzazione ottimi ha raggiunto la convergenza. Come si puo notare ii 
beamformer e in grado di annullare il segnale in corrispondenza delle direzioni d'arrivo 
degli interferenti. 

Uesame delle figg.14 e 17 mostra che applicando il metodo di equalizzazione 
oggetto della presente invenzione, la stima congiunta dei parametri ottimali di 
equalizzazione riesce a convergere verso una coppia di vettori h opt (n) e w^ r (>7) in 

grado di garantire una buona qualita nella demodulazione, anche in presenza di 
interferenti non-stazionari. Quanto appena affermato e vero fino a che Tinterferente 
non-stazionario e dello stesso ordine di grandezza del segnale utile. Se cosi non e, 
puo verificarsi un picco improvviso nel BER del segnale demodulato, difficilmente 
riducibile a valori accettabili in un tempo ragionevolmente breve rispetto alia durata 
residua del burst. 

La fig. 18 mostra tre andamenti di BER lungo il burst di altrettanti segnali 
demodulati, per un interferente che disturba fino al 50° simbolo. Le tre curve sono 

E C 
riferite a segnali d'ingresso con — r= 10 dB e tre diversi valori di — , rispettivamente 

di 0, -10, ~20dB \ gli andamenti sono mediati su 500 realizzazioni. Risulta evidente 
il peggiorare della convergenza al diminuire del rapporto C/I\ in particolare, si vede 



tar 




- 71 - 

Onlio Giu»t! 

che il BER non scende a valori accettabili fino a che il rapporto C/I resta inferiore a 
-10 dB. 

La causa di questo comportamento, e imputabile al pur basso ritardo che 
comunque affligge le decisioni a basso ritardo nella stima MLSE dei dati trasmessi. 

Proprio a causa di questo ritardo I'equalizzatore non si accorge fin da subito 
deirarrivo-di un forte interferente, di modo che, quando si rendono disponibili le prime 
decisioni provvisorie che risentono dell'interferente, idonee quindi ad indurre 
contromisure per la cancellazione dello stesso, questi avra nel contempo causato un 
gran numero di decisioni errate, di entita tale da compromettere la convergenza della 
stima effettuata dal blocco FUNZ. 

L'inconveniente viene di molto attenuato da una variante di metodo, che cerca 
di annullare nelPequalizzazione adattativa la visione ritardata dell'interferente, alio 
scopo di cancellare rapidamente forti interferenti non-stazionari. 

La variante comprende una fase di continuo monitoraggio dell'energia ricevuta 
sulle antenne, e di rivelazione di eventuali picchi superiors ad una soglia 
predeterminata. La rivelazione di un picco ad un istante 77 viene sfruttata per 
selezionare automaticamente il funzionamento dell'equalizzatore secondo le modalita 
della variante. Poiche i disturbi che devono essere rivelati sono solo quelli 
sufficientemente forti, non occorre un metodo di rivelazione dei picchi di energia 
particolarmente affinato. E possibile, ad esempio, utilizzare un metodo di rivelazione 
non coerente, noto al tecnico del ramo. Nel ricevitore di fig. 12, per semplicita, non e 
stato evidenziato un apposito blocco che svolge le suddette funzioni. 
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La fase di rivelazione e utilmente svolta in parallelo al normale funzionamento 
fuori midambolo. La variante esaurisce rapidamente i suoi effetti entro un transitorio di 
durata AT (di circa una ventina di campioni) necessario, come si vedra, 
all'assestamento nei calcoli di alcune matrici utilizzate dal blocco FUNZ. Alio scadere 
di AT conviene riselezionare automaticamente il funzionamento originale, che risulta 
piu efficiente a transitorio esaurito. 

f 

La strategia posta in atto dalla variante per poter affrontare interferenti 
stazionari di forte intensita, e quella di anticipare all'istante n la futura conoscenza 
sull'interferente da cancellare. In tal modo I'interferente e "visto" in anticipo nei calcoli, 
prima cioe che siano disponibili le prime decisioni in merito, e ralgoritmo adattativo di 
minimizzaziohe dell'errore pud cominciare ad aggiomare i coefficienti del beamformer 



per cercare di eliminare il disturbo. La suddetta strategia viene attuata nei transitorio 
AT, scomponendo in due fasi elementari (1. e 2.) Taggiornamento della matrice 
R(/?) . In pratica si agisce come segue: 

1. durante questa prima fase, che si esaurisce entro il tempo di simbolo If, viene 
anticipato I'aggiornamento della sola matrice di covarianza spaziale con i 
vettori y(w) ottenuti dai campioni dei segnali sulle antenne, presenti a partire 
dall'istante « fino all'istante IT + AT nei burst memorizzato. Viene in tal modo 
estesa la base statistica della suddetta matrice, che ora contiene sufficiente 
informazione dell'interferente da cancellare. La matrice non viene quindi piu 
toccata fino al termine del transitorio. Si passa successivamente alia fase due. 
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2. Durante questa seconda fase, per 77 < n < Ji + AT , vengono aggiornate le 
matrici di covarianza sui dati con il vettore x l {n) % e la matrice di 

cross-covarianza R xv con i vettori x'(/7) e y(//) a valle del blocco RIT. Riguardo al 
vettore x'(n) occorre precisare che a partire da 77, una volta scaduto il termine 

~~" detfe-decisiuni a basso ritarclo, edln mancanzaU i uppoitLftin^'owe^ 
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verranno tra breve precisati, questo vettore potrebbe esprimere decisioni deUtt^^ 
errate a causa del compimento della fase 1. fzs£i 
Raggiunto I'istante 77 + AT la stima delle suddette matrici procede neV0ic&o^ (/ 
convenzionale. II processo di aggiornamento della matrice R(;z) , nell'intervallo 

temporale [77, 77 + At], pud quindi essere espresso come segue: 



R( W ) = 



(77) 



dove viene esplicitato il diverso intervallo temporale utilizzato per il calcolo di R 



rispetto agli altri termini. 

L'aggiornamento "anomalo" della matrice R(/;) genera un transitorio di 
assestamento di durata AT nel calcolo delle matrici R^ e R „,. Intervallo questo che 
coincide appunto con Tampiezza dello spostamento temporale nell'aggiornamento 
della matrice di covarianza spaziale R vy . 

La stima congiunta dei parametri ottimali di equalizzazione risultera quindi 
"viziata" da questo tipo di approccio, che comunque si rivela necessario per far fronte 
a rapide variazioni delle caratteristiche di canale, consentendo al beamformer BMF di 
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introdurre comunque degli zeri in corrispondenza degli echi del nuovo interferente, 
che altrimenti non riuscirebbe affatto ad introdurre, o introdurrebbe con eccessivo 
ritardo. 

Come conseguenza del tempo di assestamento della matrice R(/;) potrebbero 



cui un'altra caratteristica della variante di metodo, e quella di aggiornare durante il 
transitorio il solo vettore w opt (n), dei pesi del beamformer BMF, mantenendo costante 
per tutta la durata del transitorio il vettore dei campioni della risposta airimpulso ad un 
valore h che esso aveva all'istante 77. La scelta premia i coefficient! del 
beamformer BMF poiche risultano meno sensibili agli errori di decisione rispetto ai 
campioni della risposta impulsiva. 

Per quanto detto, siamo ora in grade di stimare anche durante il transitorio 

[77, 77 +AT] i parametri ottimali di equalizzazione. In particolare, h opt e costante e la 
(51), ricordando la (77), consente la determinazione del vettore w 0pt , soluzione 
dell'equazione: 

(77 + AT)yv opt (n) = KXn)h opI (78) 
in cui R yy (n+AT) e calcolata all'istante 77 e poi mantenuta costante per tutto AT, 
mentre R* (77) e calcolata ad ogni istante A?-esimo del transitorio. Come si pud notare, 



a causa di h opt costante non e piu necessario risolvere le (67) e (68). 



Le precedent! considerazioni sulle modalita di aggiornamento della matrice 

R(//) della (77), nell'intervallo temporale [/F, « + AT], si trasportano pari pari 

aH'aggiornamento del fattore di Cholesky della stessa. 

Riguardo a quest'ultimo punto, ad un istante n di rivelazione del picco di 

- g n ft rq ia -r i cevut a -s ulle-ante^ = 



- viene aggiornata , entro il tempo di simbolo n f la sola matrice triangolare L con i 
vettori y(n) ottenuti dai campioni dei segnali sulle antenne, presenti nel burst 
memorizzato, a partire dall'istante 77 fino airistante ;T + AT\ 

- vengono aggiornate, per n < n < IT + AT f le matrici U(?/) e L 5 (/?) 
rispettivamente rettangolare e triangolare, con il detto vettore di aggiornamento g n 
e viene risolta I'equazione (56), mantenendo costante il vettore dei campioni della 
risposta airimpulso h opt al valore assunto al detto istante n ; 

- vengono ripristinate, alio scadere del tempo JI + AT, le modalita di calcolo 
precedent! il detto istante n . 

Per quanto detto, il fattore di Cholesky in questo intervallo temporale puo 
essere decomposto in forma triangolare come segue: 



L(») = 



L(/T + AT) 0 
U(w) L,(/i)_ 



(79) 



Con I'insieme delle considerazioni di cui sopra, relative alia possibility di 
utilizzare nei calcoli o direttamente la matrice R(«), o preferibilmente il fattore di 
Cholesky L(n) della stessa, i diversi passi della stima adattativa del solo vettore dei 
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pesi w 0 ,(a/) che caratterizzano la variante di metodo, sono ricondotti a quelli 



corrispondenti dell'algoritmo adattativo per quanto concerne la modalita di utilizzo nei 
calcoli dei dati decisi a basso ritardo {x\}. Detto algoritmo adattativo, a sua volta, 
faceva riferimento alia descrizione dettagliata dell'algoritmo svolto dal blocco FUNZ 



modifiche "relative all'avanzamento temporale AT della matrice R yy (n+AT) o del 
relativo fattore di Cholesky L(w AT) . 

Cio posto, il vettore dei pesi w opt (ti) puo essere calcolato o risolvendo la (78) o, 
preferibilmente, mediante la decomposizione in fattori di Cholesky e successiva 
risoluzione della (56) che in questo nuovo contest© diventa: 

V(JI + AT)w opt (n) = U\?i)h opt (56') 

Riguardo a questa seconda opportunita, occorre pero analizzare come 
aggiomare il fattore di Cholesky, quando esso assume I'espressione della (79), 
applicando la trasformazione unitaria Q(n) definita nella (76). Nel caso di utilizzo 
della variante, a differenza delle condizioni che valevano per le (61) e (76), risulta 
necessario aggiomare, alHstante n, il solo fattore L della matrice di covarianza 
spaziale R . Se analizziamo il meccanismo con cui si applicano le trasformazioni 
unitarie, possiamo renderci conto di come sia possibile aggiomare la sola matrice L 
"ritardando" Taggiornamento degli altri fattori. Era stato detto in precedenza, con 
riferimento alia (61) che una rotazione unitaria viene in genere attuata come una 
sequenza di trasformazioni elementari, Q 0 >Qi>'", usualmente rotazioni di Givens. II 



ritardo di aggiornamento puo essere introdotto ricordando quanto espresso nelle (62) 
e (63), dove la singola trasformazione elementare Q, agisce sulla singola colonna 
j ( /7 _l) del fattore L e sulla colonna del generatore g„. Se si applica la 
trasformazione Q, ai primi M elementi delta /-esima colonna, dove M e I'ordine d#a\ 
matrice L, e possibile ritardare I'aggiornamento dei fattori U e L Jt memorizzandp 0 1e 
trasformazioni avvenute in L e riprendendo I'aggiornamento di U e L,, unayafa^ 
che le decisioni a basso ritardo §iano ritomate affidabili, ovvero alia fine del periodo di 
^) assestamento AT di queste matrici. 

La fig. 19 mostra due andamenti di BER lungo il burst di altrettanti segnali 
demodulati, per un forte interferente che disturba fino al 50° simbolo. Un primo 
andamento (a. tratteggio) e relativo alia modalita adattativa fuori midambolo del 
metodo di equalizzazione in oggetto, mentre il secondo andamento (in solido) e 
relativo alia sua variante per forti interferenti. II segnale d'ingresso presenta 
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— = 10 dB e — = -20rfB; gli andamenti sono mediati su 500 realizzazioni. Si puo 

#0 1 



chiaramente notare che la variante di metodo inizia prima I'adattamento, riuscendo 
cosi a contenere il picco nell'andamento della probability d'errore, ed a garantire 
quindi un significativo miglioramento nella qualita del segnale demodulate. 

Ulteriori riferimenti bibliografici utili alia comprensione della "Displacement 
Theory" applicata al canale tempo-variante sono riportati nei seguenti articoli: 
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o "Time-Variant Displacement Structure and Triangular Arrays", autori: A. H. Sayed, 
H. Lev-Ari e T. Kailath, pubblicato in IEEE Trans, on Signal Processing, voL 42, No. 
5, May 1994. 

o "Fast Computation of Channel-Estimate Based Equalizers in Packet Data 

Transamission", autori: N. Al-Dhahir, J. M. Cioffi, pubblicato in IEEE Trans, on 

Signal Processing, vol. 43, No. 11, November 1995. 
o "A Computationally Efficient FIR MMSE-DFE for CCMmpaired Dispersive 

Channels", di N. Al-Dhahir, pubblicato in IEEE Trans, on Signal Processing, vol. 45, 

No. 1, January 1997. 

Gli insegnamenti del metodo di equalizzazione oggetto della presente 
invenzione, e della sua variante, trovano rispettiva applicazione negli equalizzatori 
compresi nelle schematizzazioni a blocchi funzionali delle figg.6 e 12. La 
realizzazione dei suddetti blocchi e nota al tecnico del ramo, cost come la 
realizzazione delle funzioni indicate nel testo che non trovano corrispettivi blocchi 
nelle figure, ad eccezione del blocco FUNZ e dei relativi blocchi interni PESI e CAMP. 
Inoltre, in conseguenza della diversa interrelazione tra i blocchi che una stima 
congiunta dei parametri w opt ,h opt impone rispetto ad una stima disgiunta, anche 

I'architettura deirequalizzatore non risuita nota. 

Nella descrizione dettagliata del blocco FUNZ che seguira e conveniente riferirsi 
alia sola schematizzazione di fig. 12, che consente di realizzare un unico equalizzatore 
in grado di funzionare sia durante il midambolo, che fuori, in quest'ultimo caso 
secondo la modalita adattativa, o secondo quella della variante, a seguito di 
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opportune selezioni hardware/software pilotate da eventi quali la fine del midambolo. 
o I'improvviso sopraggiungere di un forte interferente. 

Per quanto concerne la realizzazione circuitale del blocco FUNZ, e preferibile 
ricorrere a strutture circuitali note con la definizione di "array sistolici", illustrate nelle 
figg. 21, e 22. Queste strutture sono dotate di un alto grado di parallelismo, che li 
rende gia -di per se idonee ad operare sulle matrici, ma soprattutto a compiere quasi 
simultaneamente un numero di calcoli altrimenti sequenziali. Una realizzazione basata 
sull'utilizzo di un microprocessore sarebbe assai meno efficiente, proprio perche un 
comune DSP e governato da una sequenza temporale di istruzioni elementari che lo 
rendono piu idoneo a risolvere problemi seriali, a differenza dell'array sistolico che 
svolge tutti i calcoli in parallelo in un solo tempo di simbolo di 3,69 \xs (o in un tempo 
doppio se si opera congiuntamente sui semiburst). Anche volendo ricorrere ad un 
DSP di tipo RISC, questi non avrebbe un parallelismo neppure confrontabile con 
quello dell'array sistolico, che sembra quindi essere "problem oriented". Nella pratica, 
il blocco FUNZ, cosl come i restanti blocchi dell'equalizzatore di fig. 12, sono 
vantaggiosamente realizzati mediante circuiti integrati digitali di tipo ASIC (Application 
Specific Integrated Circuit). 

Con riferimento alia fig.20 si pud notare che il blocco di stima FUNZ comprende 
cinque array sistolici diversamente configurati. rispettivamente indicati con SYSTOL1, 
SYSTOL2, SYSTOL3, SYSTOL4, e SYSTOL5. Ad un primo ingresso del blocco 
SYSTOL1 perviene il vettore x(/?) riferito ai campioni della sequenza di training 
{x„}memorizzata nel ricevitore, ad un secondo ingresso perviene il vettore y(«) .' e 
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ad un terzo ingresso perviene un segnale di controllo MID che risulta vero solo in 
corrispondenza del midambolo. Su un primo bus d'uscita e presente il fattore di 
Cholesky L(n) indicato nella (61), che raggiunge un primo ingresso di un selettore a 
due ingressi SEL, mentre un secondo bus d'uscita CONT connette SYSTOL1 ad un 
primo ingresso di SYSTOL2. 

Ad un secondo ingresso del blocco SYSTOL2 perviene il vettore x'(») riferito ai 

dati della sequenza decisa a basso ritardo {x'„}; ad un terzo ingresso perviene il 
tf) vettore y(n) uscente dal blocco di ritardo RIT; un quarto ingresso BURST e un bus 

che accede in lettura ad una memoria in cui il ricevitore trascrive i campioni ></?) del 
burst di ricezione; ad un quinto ingresso perviene il coefficiente di oblio A , ad un 
sesto ingresso perviene un segnale PIC, di cui si dira; ed infine ad un settimo 
ingresso perviene il segnale di controllo MID. Sul bus d'uscita e presente il fattore di 
Cholesky Z'(n) come indicato nella (76), che raggiunge il secondo ingresso del 
selettore SEL, al cui ingresso di selezione perviene il segnale di controllo MID. II bus 
ffi uscente dal selettore SEL e connesso ad un capo di un blocco di interruzione SWT 

controllato da un segnale PIC* di durata AT proveniente da un blocco PICSEL. 
Quest'ultimo blocco riceve su tre rispettivi ingressi il segnale MID e due altri segnali 
indicati con PIC e TRANS, e sulla base di questi genera il segnale PIC. II segnale 
TRANS e attivo per un tempo AT pari alia durata del funzionamento transitorio 
previsto per la variante di metodo; come tale esso e un parametro configurabile 
dell'equalizzatore, alia stregua del coefficiente di oblio k . E sottintesa la presenza di 
una logica che provvede alia generazione ed alia gestione di tutti i segnali indicati. 
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II segnale PIC proviene da un blocco PICRIV che riceve il segnale di ricezione 
ad alta frequenza direttamente dall'array di antenne ARY. II bus uscente dall'altro 
capo del blocco SWT e connesso ad un ingresso di SYSTOL3, da cui fuoriesce il 
vettore intermedio u ;i diretto verso un primo ingresso di SYSTOL4, il cui altro 
ingresso e connesso all'ingresso di SYSTOL3. Ad una prima uscita di SYSraC 



presente un vettore h opl (n-\) diretto verso un ulteriore ingresso di SYSt^)L3. Ad 




una seconda uscita di SYSTOL4 e presente il vettore h (w) diretto verso ^^ijtrt^ 5 ^ 
ingresso di SYSTOL5, il cui secondo ingresso e connesso all'uscita del blocco SEL. 
All'uscita di SYSTOL5 e presente il vettore w^m). I due blocchi SYSTOL3 e 
SYSTOL4 sono compresi in un blocco funzionale AUTOVET, il cui ingresso coincide 
con quello di SYSTOL3 e la cui uscita coincide con quella di SYSTOL4. 

L'architettura del blocco FUNZ di fig.20 consente un pieno funzionamento 
dell'equalizzatore di fig. 12 sia durante il midambolo che al di fuori, secondo il metodo 
oggetto d'invenzione e la sua variante. 

Durante il midambolo il segnale MID risulta vero, questa condizione viene 
sfruttata per abilitare I'array SYSTOL1 e disabilitare SYSTOL2, per controllare il 
selettore SEL alia scelta di L(/?) , e determinare la chiusura deirinterruttore SWT. Di 
conseguenza il bus di SYSTOL1 e prolungato verso il blocco AUTOVET che estrae la 
parte relativa alia matrice L, e la utilizza per risolvere le (68) e (67), determinando il 
vettore h opt (n)< Piu precisamente, I'array SYSTOL3 calcola il vettore intermedio u rt 
secondo la (68) applicando il metodo della "back-substitution", mentre Tarray 
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SYSTOL4 calcola h opl (n) secondo la (67), ancora mediante "back-substitution". 
L'array SYSTOL5 estrae dal bus proveniente da SEL la parte relativa alle matrici 



L e U che utilizza, assieme al vettoreh „,,(») , per risolvere la (56) e determinare 



il 



vettore w opt («) , sempre mediante "back-substitution". 

Durante il tunzionamento adattativo, il segnale MIU neyalu disabilila laiiay 
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SYSTOL1 ed abiiita SYSTOL2, e comanda il selettore SEL alia scelta di L'(«). In 
assenza di un cambiamento di condizione logica nel segnale PIC', I'interruttore SWT 
resta chiuso ed il bus di SYSTOL2 viene prolungato verso i blocchi AUTOVET e 
SYSTOL5; questi ultimi trattano le rispettivi componenti matriciali del fattore di 
Cholesky L'(«) alia stregua di quanto detto per L(«) . 

Prima che l'array SYSTOL2 inizi ad operare, la transizione da vera a negato del 
segnale MID comanda il trasferimento del contenuto dell'array SYSTOL1 sul bus 
CONT per essere ricopiato entro le celle di SYSTOL2. La suddetta operazione 
consente I'aggiomamento del preesistente fattore di Cholesky L(«) secondo la nuova 
^ modalita di calcolo che prevede I'utilizzo del coefficiente di oblio X. Se non fosse stato 

previsto il trasferimento di cui sopra, l'array SYSTOL2 stimerebbe ex-novo il fattore 
L'Cw) ma con minori possibility rispetto alia situazione precedente. 

Quando il blocco PICRIV rivela, ad un istante 77 , la presenza di un picco di 
energia sulle antenne superiore alia soglia prefissata, il blocco PICSEL discrimina tra 
due possibili condizioni di funzionamento. Una prima condizione e quella di forte 
interferente in presenza di midambolo, in tal caso il segnale PIC risulta falso e la 
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stima durante il midambolo non viene disturbata. La seconda condizione e quella di 
forte interferente al di fuori del midambolo, in tal caso il segnale PIC vero comanda 
I'apertura deirinterruttore SWT, e controlla opportuni commutatori inclusi nell'array 
SYSTOL2 che apportano una modificazione provvisoria nella struttura dell'array tale 
da renderlo idoneo ad operare come previsto nella variante di metodo. Con PIC vero 




j| blocco A'UTOVET non riceve piu nulla in ingresso e provvede a congelare il valore 
del vettore h opt all'istante n fino alio scadere di AT . II blocco SYSTOL5 risolve la 
(56), sempre mediante "back-substitution", utilizzando il vettore h opt e le matrici 
L e U del fattore di Cholesky £•(>/) che deriva dalle suddette modifiche strutturali 
apportate all'array SYSTOL2, determinando in tal modo il vettore w opt (>?) 

Alio scadere del tempo 77 + AT la nuova condizione del segnale PIC comanda 
la chiusura deirinterruttore SWT, ma la configurazione dell'array SYSTOL2 rimane 
fino alia fine del burst quella modificata, ed il funzionamento dell'equalizzatore ritorna 
a quello del metodo adattativo. 
^p) Per quanto concerne la realizzazione circuitale dei diversi blocchi di fig.20, 

occorre precisare quanto segue: 

- i blocchi PICRIV, PICSEL, RIT, SEL, e SWT, sulla base della descrizione 
funzionale fornita, sono realizzabili applicando le comuni conoscenze tecniche; 

- gli array sistolici SYSTOL1, SYSTOL3, SYSTOL4, e SYSTOL5 sono 
schematizzabili da un'architettura di celle comune a quella di tutti gli array sistolici 
noti di forma triangolare. La loro descrizione richiede pertanto di specificare la 
natura delle celle, per evidenziare dove gli array si diversificano tra loro e da quelli 



dell'arte nota. Una descrizione di tipo funzionale delle celle che indichi il tipo di 
operazioni aritmetiche compiute sul segnale entrante e sul contenuto preesistente 
per ottenere il nuovo contenuto ed il segnale uscente, e da ritenersi sufficiente alia 
loro realizzazione sfruttando le comuni conoscenze tecniche nel campo degli array 
sistolici e dei gate array. 
I I'array sistolico SYSTOL2 si avvale della medesima struttura di SYSTOL1, ma 
include delle memorie FIFO e dei commutatori elettronici nella catena di discesa, 
come si vedra in fig.22. Le modalita di ottenimento in chip delle suddette FIFO e 
dei relativi commutatori, sono ancora quelle che consentono di ottenere i diversi 
array, vale a dire la tecnologia dei gate array di tipo ASIC. 

Ovviamente, il blocco FUNZ pud essere realizzato mediante architetture 
differenti rispetto a quella indicata, e possibile ad esempio utilizzare un unico array 
sistolico che conglobi le funzioni di SYSTOLI e SYSTOL2, intal caso risulterebbe 
superfluo il trasferimento del contenuto alia fine del midambolo. Questo nuovo array 
deve comunque poter modificare alia fine del midambolo le operazioni aritmetiche che 
definiscono il contenuto delle celle, in quanto occorre ora tenere conto del coefficiente 
di oblio k . 

Con riferimento alia fig. 21, in cui I'array sistolico SYSTOL1 viene schematizzato 
mediante blocchi funzionali, si pud notare che ci sono singoli blocchi a forma quadrata 
tra loro interconnessi, indicati con l i , JI i e l s . , e blocchi periferici di forma circolare 
c i connessi ai primi. La struttura dell'array e di forma triangolare, cosi come quella 
della matrice che esso riproduce. In effetti ciascun blocco I , il, , l s . , e 
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rappresenta un elemento cli questa matrice. Come si pud Hgtale, i : 
raaaruppati nella zona del vertice acuto superiore corrisporiqfente ^ll|i'.posii»bne della 
matrice triangolare L, i blocchi l t _.,c, sono raggruppati nella zona del vertice acuto 
inferiore corrispondente alia posizione della matrice triangolare L, , ed infine i blocchi 
jj- ... W n aa r, r pnt i pniin rimnnnntft zona di fnrma ""adrata f n p j u oenera l mente . 
rettangolare) corrispondente alia posizione della matrice U . I blocchi circolari c. 
costituiscono la diagonale di chiusura dell'array, ed i due blocchi alle estremita 
indicano I'inizio e la fine dell'aggiornamento. Per comodita di rappresentazione sono 
indicate matrici di 3x3. 

La fig.22 indica I'architettura dell'array SYSTOL2, ottenuta da quella di fig.21 
interrompendo le linee d'uscita discendenti dei blocchi corrispondenti agli elementi 
dell'ultima riga della matrice triangolare L , inserendo delle identiche memorie FIFO di 
lunghezza tale da introduce un ritardo pari al transitorio AT , e connettendo le uscite 
delle FIFO a rispettivi blocchi corrispondenti agli elementi della prima riga della 
matrice U . L'array SYSTOL2 e in grado di funzionare secondo la modalita adattativa 
e secondo la variante di metodo, verranno quindi illustrati i mezzi che consentono 
questa possibilita. 

La fig.23 mostra un commutatore elettronico COM, identico ad altri inclusi 
nell'array SYSTOL2, inserito tra due punti indicati con A e B in fig.22, corrispondenti 
alle porte d'ingresso e uscita di una generica FIFO. II punto centrale del commutatore 
COM e connesso al punto A di provenienza del segnale da commutare, mentre la 
FIFO e connessa tra i due terminali d'uscita di COM, uno dei quali e connesso al 
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I punto B. II segnale PIC attivo comanda la commutazione in modo da inserire la FIFO 

tra A e B come mostra la posizione di fig. 23; viene in tal modo ritardato di AT il 
segnale che si propaga sulla corrispondente colonna. Alia fine del transitorio AT, 
tuttavia, la FIFO rimane inserita fino alia fine del burst, istante in cui avviene la 
commutazione opposta che cortocircuita i punti A e B, ed il ritardo viene escluso 
ritornandq al funzionamento adattativo in presenza di un nuovo burst. In figura viene 
*|j pure evidenziato un segnale di orologio CKT a frequenza di simbolo che temporizza il 

funzionamento della FIFO. Ci saranno tanti blocchi COM quante sono le colonne 
sovrastanti le FIFO. 

Con riferimento alia fig.24, si pud notare che i blocchi circolari sono interessati 
da un segnale entrante in e da un segnale uscente p. Anche i blocchi quadrati sono 
interessati dai medesimi segnali, ma jn aggiunta hanno un ulteriore segnale entrante 
pe un ulteriore segnale uscente out. Fanno eccezione i blocchi quadrati deirultima fila 
orizzontale e I'ultimo blocco circolare della stessa, ove manca il segnale uscente p. II 
contenuto dei blocchi, identicamente per ciascun tipo di blocco, viene indicato con a. 

Gli array SYSTOL1 e SYSTOL2 sono disposti nelle figure in modo che i segnali 
entrant!, corrispondenti al vettore di aggiornamento: 



%n = 



x(/i) 



entrino da sinistra. Piu precisamente, i segnali numerici y x (n\...,y£?i\....,y M (n) 
sono diretti verso i blocchi della prima linea verticale dell'array, corrispondenti agli 
element! della prima colonna della matrice triangolare L , mentre i segnali numerici 
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x, (/?),..., x i (n\....,x L (?i) sono diretti verso i blocchi corrispondenti agli elementi della 
prima colonna della matrice (rettangolare) U . Per quanto concerne I'array SYSTOL2, 
il vettore x(n) sia durante il funzionamento adattativo che durante la variante coincide 
con il vettore x'(n)\ il vettore y(//) durante il funzionamento adattativo sara quello 



y(n + AT) costituito dai campioni y } (n + AT), ,(/7 + AT\.. ..,y M (ti + AT) del 
burst memorizzato, trasferiti in SYSTOL2 ed ivi utilizzati nell'aggiornamento anticipato 
del fattore di Cholesky L(;7 + AT) della matrice di covarianza spaziale R >y (ii + AT) . 

Anche se non esplicitamente indicato, e previsto un segnale di orologio 
generale dell'array delle figg.21 e 22, a frequenza molto piu alta di quella di simbolo, 
che perviene simultaneamente a tutti i blocchi l t , ti i , l s . , e c i per controllarne il 
singolo funzionamento; il banco di memorie FIFO di fig.22 viene aggiornato a tempo 
di simbolo. 

La disposizione geometrica dei blocchi illustrata nelie figg.21 e 22 viene 
rispettata sull'area dei chip, ove i blocchi costituiscono altrettante celle; cio consente 
al segnale che origina nel vertice superiore di propagarsi lungo la catena discendente. 
E evidente che il corretto contenuto delle celle sara disponibile solo dopo un tempo 
sufficiente alia completa propagazione del segnale verso restremo opposto della 
diagonale costituita dai blocchi circolari, seguendo il percorso piu lungo tra le celle. I 
singoli segnali entranti nelle parte sinistra dell'array sono associati ad un 
corrispondente numero di frecce che indica il ritardo con cui si esplica il loro effetto. 
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Nel funzionamento, gli array sistolici SYSTOL1 e SYSTOL2 sono 
rispettivamente riferiti all'aggiornamento del fattore di Cholesky Z(n) della matrice 
R(/?) come indicato nella (61) e nella (76), che per comodita riproponiamo: 

[L(») o] = [L(//-l) g„]Q(») (61) 
[L(A/rT]~[^L(7i 1) fe n ]uO') ^ 




Durante i forti interferenti I'aggiornamento indicato dalla (76) viene modificato 
dall'introduzione delle FIFO come visto nelle figg.22 e 23, ma le relazioni che 
verranno ora indicate continueranno a valere sugli elementi della matrice L(w + AT) . 

La rotazione unitaria Q(;?) viene attuata associando a ciascuna cella una 
relativa operazione compiuta sui segnali entranti e sul contenuto preesistente, per 
ottenere i segnali uscenti ed il nuovo contenuto. Nella fattispecie, ad esempio, della 
(76): 

- i segnali in si propagano in senso orizzontale e corrispondono ai segnali 
y^nl- ^M >yM(») e x l (n),...,x i (nl ...,x L (") all'ingresso delle celle della 
prima colonna dell'array; per le restanti celle indicate dai blocchi quadrati i segnali 
in corrispondono ai segnali out uscenti dalle celle precedenti nel senso di 
propagazione orizzontale del segnale; 

- i segnali p si propagano in senso verticale e vengono generati all'uscita di 
ciascuna cella, sia quadrata che circolare, calcolando la seguente espressione: 



(80) 



ccJA 



dove X e il coefficiente di oblio, e a e il nuovo contenuto delle celle calcolato 
come segue, relativamente alle sole celle circolari corrispondenti alia diagonale 
principale della matrice L : 



a <- yjx\a\ 2 +|///| 



- j l-nvjo vo-flf-Q ssond o-p an ta to~d a 1 1 tH re nci at 




per le restanti celle quadrate il nuovo a viene determinate OQ^^tamerrte al 
segnale out che si propaga in senso orizzontale, vale infatti la seguente relazione 
matriciale: 



[a out\ [a^fx /w] 



1 p" 

p -1 



(82) 



dove il vettore puntato dalla freccia indica i nuovi valori di a e out . 

Relazioni analoghe alle (80), (81), e (82) valgono per I'aggiornamento secondo 
la (61), ponendo X = 1. 

Per quanto concerne la fig.22, I'introduzione delle memorie FIFO nella catena di 
discesa dell'array sistolico e la semplice traduzione circuitale di quanto detto 
illustrando la variante di metodo a proposito dell'aggiornamento del solo fattore di 
Cholesky L(77 + AT) . Si e infatti detto che e possibile ritardare I'aggiornamento dei 
fattori U e L, memorizzando i risultati delle trasformazioni elementari Q.. II ritardo 
introdotto all'istante n viene poi recuperato all'istante n + AT t cioe quando lo 
stimatore STIM e in grado di fornire i simboli decisi affidabili e pub quindi riprendere 
un aggiornamento normale del fattore di Cholesky L(/?) . 
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Per quanto concerne la descrizione funzionale delle celle appartenenti agli 
array sistolici SYSTOL3, SYSTOL4 e SYSTOL5, ovvero le relazioni matematiche 
sviluppate nelle celle sulla base del contenuto precedente e dei nuovi ingressi per 
generare delle nuove uscite, il tecnico del ramo e posto in grado di ottenere i suddetti 
array sistolici, applicando i passi del metodo di back-substitution descritto nel testo per 



risolvere le relative equazioni indicate, e sfruttando le comuni conoscenze tecniche 
nel campo degli array sistolici. 

Vantaggiosamente, I'equaiizzatore appartenente al ricevitore di fig. 12 (in pratica 
I'intera struttura a blocchi mostrata in figura) viene duplicato onde poter processare 
parallelamente i due semiburst ai due lati del midambolo. Un circuito di ricostruzione 
del burst posto a valle dei due equalizzatori pone neiresatto ordine temporale i dati 

delle due sottosequenze stimate, ottenendo un'unica sequenza (x n }dei dati 

trasmessi stimata a massima verosimiglianza. 



m 
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RIVENDICAZIONI 

1. Metodo di equalizzazione per la cancellazione di interferenti 
isofrequenziali in segnali ricevuti da antenne a schiera (ARY) di stazioni radio base di 
un sistema radiomobile cellulare, detti segnali (y(/?)) essendo ricevuti sotto forma di 
burst di simboli di un burst di trasmissione, comprendente una sequenza di t raining 




nota al ricevitore . detti burst ricevuti essendo demodulati, digitalizzati, 
memorizzati, e sincronizzati mediante discriminazione della sequenza di training, 
caratterizzato dal fatto che comprende le seguenti fasi: 

a) filtraggio spaziale dei segnali provenienti dalle antenne (y(//)) , utilizzando un filtro 

spaziale (BMF) che combina tra loro sulla base di opportuni pesi {w opt (n)) dei 
simboli corrispondenti nei diversi burst memorizzati; 

b) filtraggio temporale del segnale (y out (n)) generato da detto filtro spaziale (BMF) 

utilizzando un filtro adattato (FADAT) a campioni (h opt (n)) della risposta impulsiva 
del canale trasmissivo; 

c) ottimizzazione congiunta di detti pesi e campioni (w 0/> , (/?), h 0/J ,(/7)) sulla base dei 

simboli di detta sequenza di training ({*„}) e dei detti segnali provenienti dalle 
antenne (y(n)) presenti nei burst memorizzati; 

d) stima dei simboli ({*„}) del burst di trasmissione corrispondenti alia sequenza 
uscente dal filtro temporale adattato (FADAT); 



e) aggiornamento dei valori ottimizzati di detti pesi e campioni [w opt (//), h opt (//)), al di 



fuori della finestra temporale della sequenza di training, mediante ripetizione delle 
precedent! fasi a)...d), sostituendo nella fase c) i simboli della sequenza di training 



equa li zzaz i one adattativa. = 
2. Metodo di equalizzazione secondo la rivendicazione 1, caratterizzato dal 

fatto che detta ottimizzazione congiunta di pesi e campioni (w opt (ti), h o/7f (>0) 6 



compiuta minimizzando un funzionale (FUNZ) corrispondente all'errore quadratico 
medio che nasce dal confronto tra corrispondenti campioni /7-esimi di due sequenze, 

di cui, una prima e ottenuta sottoponendo i simboli della sequenza di training ({*„}) 
ad un filtro temporale (CAMP) che stima i detti campioni della risposta impulsiva del 
canale trasmissivo (h(/?)) ( ed una seconda sequenza e ottenuta sottoponendo i 
campioni /7-esimi dei detti burst memorizzati {yin)) ad un filtro spaziale (PESI) che 

stima i detti pesi (w(/?)) . 

3. Metodo di equalizzazione secondo la rivendicazione 1 o 2, caratterizzato 
dal fatto che la detta ottimizzazione congiunta di pesi e campioni converge entro una 

finestra temporale della detta sequenza di training ({*„}) di lunghezza inferiore o 

uguale a 26 intervalli di simbolo. 

4. Metodo di equalizzazione secondo una qualunque delle rivendicazioni 

precedent!, caratterizzato dal fatto che quando la detta sequenza di training ({*„}) e 




rendendo con cio detta 




/ 



inclusa in un midambolo del burst di trasmissione, la detta equaliziazione inizia^alle 
estremita del midambolo e procede parallelamente nei due sertsti^^bdrst, 
preventivamente memorizzato e sincronizzato. 

5. Metodo di equalizzazione secondo una qualunque delle rivendicazioni da 

1 a 3, caratterizzato dal fatto che quando la detta sequenza di training ({*„}) e 

inclusa in un midambolo del burst di trasmissione, la detta equalizzazione inizia da 
una prima estremita del midambolo e procede verso la fine del burst, preventivamente 
memorizzato e sincronizzato, e successivamente riprende dalla seconda estremita del 
midambolo verso I'altra estremita del burst, usando i dati stimati del primo semiburst 
per avere un insieme di dati stimati piu lungo per ('equalizzazione del secondo 
semiburst 

6. Metodo di equalizzazione secondo una qualunque delle rivendicazioni da 

2 a 5, caratterizzato dal fatto che detta ottimizzazione congiunta di pesi e campioni 
compiuta nella detta fase c) comprende i seguenti passi: 

d1) assegnazione di valori arbitrari, preferibilmente nulli, agli elementi di una matrice 



triangolare strutturata L(/z) , definita come segue: L(/?) = ^ ^ 



, dove 



L(w) e L s (n) sono matrici triangolari ed U(//) e una matrice quadrata, o piu 

generalmente rettangolare; 
d2) assegnazione di valori arbitrari, di cui almeno uno non nullo, agli elementi di un 
vettore h opt (n) che memorizza i detti campioni ottimizzati della risposta 
impulsiva del canale trasmissivo; 
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d3) assegnazione di valori arbitrari, di cui almeno uno non nullo, agli elementi di un 
vettore w opl (n) che memorizza i detti pesi ottimizzati del filtro spaziale (BMF); 



d4) aggiornamento di detta matrice strutturata L(/?) applicando all'istante n -esimo 
la seguente relazione iterativa: [E(>?) o]' = [L(«-l) g„]Q(«)> dove g„ « 11 



® 



seguente vettore di aggiornamento g n - 



strutturato in due vettori 



component!, di cui: yfc) = [j/ l (w), - K 2 00.----.-Vw(")] r e un P rimo vettore i cui M 
elementi sono simboli provenienti da M antenne all'istante n -esimo, present! in 

altrettanti detti burst memorizzati, ed x(n) = [x n ( n \ x n+\( n \ ^ x n+L-\( n )] T ® un 
secondo vettore i cui L elementi sono dei simboli memorizzati della sequenza 
trasmessa ({xj, {*'„}) . rispettivamente di training, o stimata, assumendo una 
durata della risposta impulsiva pari a L simboli; essendoQ(/7) una 
trasformazione matriciale di rotazione unitaria, attuata tramite una sequenza di 
trasformazioni elementari Q 0 , <},,•••, che applicata a g„] produce 

I'annullamento della colonna g„ introducendo uno zero per ogni trasformazione 
elementare; 

d5) determinazione di un vettore di pesi h n mediante la seguente espressione: 
L,(/?)L* (w)h« = h ^r(" _1 )^ da cui viene aggiornato come segue il detto 



vettore dei campioni ottimi h opt {?7) = . ; 
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cl6) aggiornamento del detto vettore dei pesi ottimi mediante la seguente 
espressione: V(n)w opt (") = U*(») h op , (>0 • 

1. Metodo di equalizzazione secondo la rivendicazione 6, caratterizzato dal 
fatto che detta espressione calcolata al passo d5) viene scomposta in due fasi 
elementari di calcolo, introducendo un vettore intermedio incognito u„ definito dalla 




seguente prima equazione: L*(//)h„ = u„, essendo u„ una soluzione della 
seguente seconda equaziope: L,(/?)u„ = h opt (n-\) che sostituita a sua volta 
nell'equazione di partenza fornisce h rtI semplificandone il calcolo. 

8. Metodo di equalizzazione secondo la rivendicazione 6 o 7, caratterizzato 
dal fatto che nella detta fase e) la detta matrice L(/?- l) viene moltiplicata per VX, 
dove X e una costante reale arbitraria, inferiore all'unita, scelta in relazione alia non 
stazionarieta del canale rendendo piu efficace la detta equalizzazione adattativa. 

9. Metodo di equalizzazione secondo una qualunque delle rivendicazioni 
precedents caratterizzato dal fatto che detta stima dei simboli del burst di trasmissione 

© compiuta nella fase d) e una stima di sequenza ({x n }) a massima verosimiglianza 

effettuata mediante I'algoritmo di Viterbi, utilizzando per questo i detti campioni 

ottimizzati della risposta impulsiva del canale (h opt (n)) per calcolare i valori di 

opportune metriche di ramo di un trellis che indica tutte le possibili transizioni di stato 
del ricevitore, Tespressione delle metriche di ramo essendo nota nei detti ricevitori. 

10. Metodo di equalizzazione secondo le rivendicazioni 1 e 9, caratterizzato 

dal fatto che nella fase di aggiornamento e) i simboli della sequenza di training ({*„}) 
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vengono sostituiti con una sequenza di simboli ({*'„}) stimata a basso ritardo 
forzando I'algoritmo di Viterbi (STIM) ad emettere delie decisioni alio scadere di un 
tempo inferiore a quello impiegato per decidere la detta sequenza stimata a massima 

verosimiglianza ({*„})• 
11 Metodo-di-QquafeTa-zi^^ 




fatto che, in corrispondenza di ogni nuovo simbolo della detta sequenza stimata a 

basso ritardo ({*'„}), gli £;elementi del detto vettore x(/?) sono rimpiazzati con 

altrettanti simboli ricopiati in blocco dal percorso selezionato nel trellis al tempo di 
simbolo attuale, prevenendo con cio la propagazione degli errori nella stima 

congiunta dei detti pesi e campioni ottimi (w opt (//), h opt . 

12. Metodo di equalizzazione secondo la rivendicazione 10 o 11, 
caratterizzato dal fatto che il detto vettore y(//) subisce un ritardo uguale a quello 

complessivamente subito dalla detta sequenza stimata a basso ritardo ({*'„})• 

13. Metodo di equalizzazione secondo una qualunque delle rivendicazioni 
precedenti, caratterizzato dal fatto che comprende inoltre una fase di continuo 
monitoraggio dell'energia ricevuta sulle antenne (ARY) e di rivelazione di eventuali 
picchi superiori ad una soglia predeterminata. 

14. Metodo di equalizzazione adattativa secondo le rivendicazioni 6 e 13, 
caratterizzato dal fatto che, ad un istante n di rivelazione di un detto picco di energia 
ricevuta sulle antenne vengono compiute le seguenti azioni: 



# 
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aggiornamento anticipato, entro il tempo di simbolo /? j|Selfa so: 5 * 
triangolare L della detta matrice L(//) utilizzando i vettor^j//)^ 
campioni dei segnali sulle antenne, presenti nei burst memorizzati, a partire 
dall'istante Ji fino ad un istante Ji + AT che delimita un transitorio di assestamento 
dei valori calcolati neiraggiornamento della detta matrice L(/?) ; 

- aggiornamento, per Ji < n < Ji + AT , delle matrici U(/?) e L 5 (/7) della detta. 
matrice triangolare L(/7) ^on il detto vettore aggiornamento g n ed esecuzione del 
detto passo d6) che aggiorna il detto vettore dei pesi ottimi (w op ,) , mantenendo 

pero costanti i detti campioni della risposta impulsiva del canale [h opt ) ai valori 

assunti al detto istante n \ 

- ripristino, alio scadere del tempo Ji + AT, delle modalita di aggiornamento della 
detta matrice L(??) precedent! il detto istante Ji . 

15. Metodo di equalizzazione secondo la rivendicazione 14, caratterizzato 
dal fatto che detto aggiornamento anticipato della sola matrice triangolare L , viene 
effettuato applicando la generica detta trasformazione elementare Q. ai primi M 

element! della /-esima colonna della detta matrice triangolare L(/7), doveM e I'ordine 
della matrice L, ripetendo la trasformazione elementare Q. per tutti i valori 
dell'indice / e memorizzando le trasformazioni avvenute. 

16. Equalizzatore in ricevitori di stazioni radio base di un sistema radiomobile 
cellulare, in particolare per la cancellazione di interferenti isofrequenziali in segnali 
ricevuti da antenne a schiera (ARY) di stazioni radio base di un sistema radiomobile 
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cellulare, detti segnali (y(/?)) essendo ricevuti sotto forma di burst di simboli di un 
burst di trasmissione, comprendente una sequenza di training ({*„}) nota al 

ricevitore, detti burst ricevuti essendo demodulati, digitalizzati, memorizzati, e 
sincronizzati mediante discriminazione della sequenza di training, caratterizzato dal 
-fatto~che co mprend e: 




a) un filtro spaziale (BMF) dei segnali provenienti dalle antenne (yO?)) che combina 
tra loro simboli corrispondenti dei diversi burst memorizzati sulla base di opportuni 

pesi (w^, («)) ; 

b) un filtro temporale (FADAT) adattato ai campioni (h opl (nj) della risposta impulsiva 

del canale trasmissivo, che filtra il segnale (y ou , («)) proveniente dal detto filtro 
spaziale (BMF); 

c) mezzi di ottimizzazione congiunta (SYSTOL1, 3,4,5) di detti pesi e campioni 
{w opt (n),h opl (n)) sulla base dei simboli di detta sequenza di training ({*„}) e dei 
detti segnali provenienti dalle antenne (y(«)) presenti nei burst memorizzati; 

d) un organo decisore (STIM) che stima i simboli ({*„}) del burst di trasmissione 

corrispondenti alia sequenza uscente dal filtro temporale adattato; 

e) mezzi di aggiornamento (SYSTOL2,3,4,5) dei valori ottimizzati di detti pesi e 

campioni {w opt (n), h opl (nj) che agiscono al di fuori della finestra temporale della 
sequenza di training ({*„}) e che comprendono i precedenti mezzi a)...d) in cui i 
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mezzi di ottimizzazione c) agiscono ora sulla base di simboli stimati ({* ! „}) 

dall'organo decisore d), rendendo con cio detta equalizzazione adattativa. 

17. Equalizzatore secondo la rivendicazione 16, caratterizzato dal fatto che 
detti mezzi di ottimizzazione congiunta (SYSTOL1 ,3,4,5) comprendono un primo 

come in un vettore algebrico y(/?) di M elementi, tanti quante sono le antenne, ed i 
simboli della detta sequerfza di training ({*„}) organizzati come in un vettore 




algebrico x(n) di L elementi, pari alia lunghezza della risposta impulsiva del canale, e 
genera in corrispondenza di ciascun tempo di simbolo //-esimo un insieme di valori 
contenuto in una disposizione spaziale di celle che riproduce la forma di una matrice 
triangolare strutturata L(n) i cui elementi hanno i valori dei corrispondenti contenuti 
delle celle; opportuni sottoinsiemi di detti valori presi al tempo di simbolo w-esimo, 
aggiornando la stima di un vettore h opt (??) che memorizza i detti campioni ottimi della 
risposta impulsiva e la stima di un vettore w c/Jf («) che memorizza i detti pesi ottimi. 

18. Equalizzatore secondo la rivendicazione 17, caratterizzato dal fatto che i 
detti vettori y(n) e x(n) pervengono incolonnati, nell'ordine indicato, all'ingresso delle 
celle di detto primo array sistolico (SYSTOL1) corrispondenti agli elementi della prima 
colonna di detta matrice L(w), I'insieme di celle della matrice triangolare includendo 
mezzi di calcolo della segue relazione matriciale: 
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dove: 

- a puntato dalla freccia indica il nuovo contenuto deila cella; 

- out e un segnale uscente; 

- in e un segnale entrante in ciascuna cella; e 

- p e un segnale sia entrante che uscente, calcolato come segue: 




w 

le celle corrispondenti agli e'lementi della diagonale principale di detta matrice L(;z) 
contenendo mezzi di calcolo alternative tali per cui per cui il segnale p e solo 
uscente, il segnale out e assente, ed il contenuto a e aggiornato come segue: 



a <r- tJ\o\ 2 +|/>/| 2 . 

19. Equalizzatore secondo le rivendicazioni 16 e 18, caratterizzato dal fatto 
che detti mezzi di aggiornamento (SYSTOL2 ) 3,4,5) comprendono un secondo array 
sistolico (SYSTOL2) le cui celle includono mezzi di calcolo che differiscono da quelli di 
detto primo array sistolico (SYSTOL1) per il solo fatto che il contenuto a viene 

moltiplicato per un valore numerico costante VX, dove X e una costante reale 
arbitraria, inferiore all'unita, scelta in relazione alia non stazionarieta del canale per 
rendere piu efficace la detta equalizzazione adattativa. 

20. Equalizzatore secondo la rivendicazione 18 o 19, caratterizzato dal fatto 
che detta matrice triangolare immagine di detti primo e secondo array sistolici 
(SYSTOL1 , SYSTOL2), e definita come segue: 
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in cui gli elementi L(n) e L,(w) sono pure matrici triangolari, e VOW^^fJ}9O0i t o 



quadrata. 

21 . Equalizzatore secondo la rivendicazione 20, caratterizzato dal fatto che i 
detti mezzi di ottimizzazione (SYSTOLi ,3,4,5), o i detti mezzi di aggiornamento 
(SYSTOL2,3,4,5), comprendono inoltre un terzo array sistolico (AUTOVET) ed un 
quarto (SYSTOL5) posti in cascata, detto terzo array ricevendo il detto elemento 
L, s (n) ed aggiornando gli elementi del detto vettore h opt (n) , e detto quarto array 

sistolico ricevendo i detti elementi L(/7) e U(/?) ed il vettore h opt (/?) , ed aggiornando 

gli elementi di detto vettore w opt («) . 

22. Equalizzatore secondo una qualunque delle rivendicazioni da 16 a 21, 
caratterizzato dal fatto che detto organo decisore (STIM) e uno stimatore di Viterbi che 

stima a massima verosimiglianza la sequenza di simboli ({*„}) del burst di 
trasmissione, utilizzando per questo i detti campioni ottimizzati della risposta 
impulsiva del canale h opt (n) per calcolare i valori di opportune metriche di ramo di un 

trellis che indica tutte le possibili transizioni di stato del ricevitore, Tespressione delle 
metriche di ramo essendo nota nei detti ricevitori. 

23. Equalizzatore secondo le rivendicazioni 16 e 22, caratterizzato dal fatto 
che inoltre comprende: 

- mezzi di memorizzazione (DECTEMP) dei simboli appartenenti ad un percorso di 



stima selezionato nel trellis, con un ritardo di durata AT inferiore al ritardo di 
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emissione dei simboli della detta sequenza stimata a massima verosimiglianza 
\ {{*»}) • ■ simboli memorizzati essendo organizzati come un vettore algebrico x'(n) 
che perviene all'ingresso di detto secondo array sistolico (SYSTOL2); 
- mezzi di ritardo (RIT) che ritardano di AT il detto vettore algebrico y(/?) prima che 
r*nn'""na4l-dettQ-secQndQ^ s 



- mezzi logici (MID) di rilevazione della coincidenza, o non, tra la posizione nel burst 
dei simboli ricevuti e la finestra temporale della sequenza di training nei simboli 

trasmessi ({*„}); 

- mezzi di selezione (SEL), controllati dai detti mezzi logici (MID) di rilevazione, che, 
in caso di coincidenza con la sequenza di training, prolungano il bus su cui viaggia 
il detto insieme di valori L(;?) all'ingresso di detto terzo array sistolico (AUTOVET) 
e, parallelamente, all'ingresso di detto quarto array sistolico (SYSTOL5), mentre in 
caso di non-coincidenza prolungano invece il bus che trasporta detto insieme di 

valori L,'(n) ; 

- mezzi di trascrizione (CONT) del contenuto di detto primo array sistolico 
(SYSTOL1) nelle corrispondenti celle del detto secondo array sistolico (SYSTOL2), 
comandati dai detti mezzi logici (MID) di rilevazione della sequenza di training, che 
attivano la trascrizione in concomitanza con la transizione di chiusura della finestra 
temporale della sequenza di training. 

24. Equalizzatore secondo la rivendicazione 23, caratterizzato dai fatto che 
inoltre comprende: 
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mezzi di rilevazione di un istante 77 di accadimento di un picco di energia nei detti 
segnali ricevuti sulle antenne, e di generazione di un segnale di presenza di forte 
interferente (PIC), attivo esternamente alia finestra temporale della sequenza di 
training per la detta durata A7; 

mezzi di interruzione (SWT) di detto prolungamento verso I'ingresso di detto terzo 



array sistolico (AUTOVET), controllati da detto segnale di presenza di forte 
interferente (PIC); 

- mezzi di congelamento del contenuto (h opt ("i) di detto quarto array sistolico 

(SYSTOL5) al detto istante 77, controllati da detto segnale di presenza di forte 
interferente (PIC); 

- mezzi (COM, FIFO) di variazione della configurazione interna dj detto secondo 
array sistolico (SYSTOL2), controllati da detto segnale di presenza di forte 
interferente (PIC). 

25. Equalizzatore secondo la rivendicazioni 24, caratterizzato dal fatto che i 
detti mezzi (COM, FIFO) di variazione della configurazione interna di detto secondo 
array sistolico (SYSTOL2) includono dei commutatori elettronici (COM) sezionanti le 
linee dell'array che connettono le celle corrispondenti agli elementi dell'ultima riga 
della detta matrice triangolare L(«) con gli elementi della prima riga della detta 

\ matrice, rettangolare U(/?) ? e che connettono altrettante memorie FIFO tra le due 

estremita delle linee sezionate. 

26. Equalizzatore secondo una qualunque delle rivendicazioni da 17 a 25, 
caratterizzato dal fatto che i detti primo (SYSTOL1) e secondo (SYSTOL2) array 
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sistolici sono conglobati in un unico array avente la comune struttura triangolare cli 
celle dei detti primo e secondo array, e che include mezzi aggiuntivi di moltiplicazione 
di detto contenuto a delle singole celle per il detto valore numerico costante VX , 
detti mezzi di moltiplicazione essendo controllati dai detti mezzi logici di rilevazione 
(MID) della finestra temporale della sequenza di training. 

27. Equalizzatore secondo la rivendicazione 21, caratterizzato dal fatto che 
detto terzo array sistolico (AUTOVET), e ripartito nella cascata di due array sistolici 
(SYSTOL3, SYSTOL4) di identica struttura triangolare, di cui, I'array a monte 
(SYSTOL3) riceve il detto insieme di valori L,(/?) , ed il precedente vettore h opt (n- l) 
dall'array posto a valle, e calcola un vettore intermedio u„ che fornisce all'array a 
valle (SYSTOL4), a cui perviene inoltre un insieme di valori L*(w), per il calcolo di 
detto vettore h op( (n) . 

28. Equalizzatore secondo una qualunque delle rivendicazioni da 16 a 27, 
caratterizzato dal fatto che quando detta sequenza di training appartiene ad un 
midambolo corto rispetto alia lunghezza dell'intero burst, il detto equalizzatore viene 
duplicato, i due equalizzatori processando parallelamente I'uno, o I'altro, dei due 
semiburst ai due lati del midambolo, e fornendo una rispettiva sottosequenza stimata 
ad un circuito di ricostruzione del burst che pone neiresatto ordine temporale i simboli 
delle due sottosequenze stimate, ottenendo un'unica detta sequenza dei dati 

trasmessi stimata a massima verosimiglianza {x n } . 
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